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i In this thesis we show that the effective Lagrangian models, encoding the rele- 

I vant symmetries of the underlying fundamental gauge theory for strong inter- 

im ' actions (QCD), provide a reasonable understanding of the interactions among 

light mesons at intermediate energies as well as of the properties of Heavy 
baryons. In the first part of the thesis we show that it is possible to build 
a chiral and crossing symmetric model for light meson interactions at inter- 
mediate energies, which we have called the Chiral Resonance Model (ChRM). 
. Using the previous scheme we can understand tttt— scattering up to the 1 GeV 

region by considering the resonance exchange together with the contact term 
contributions. We also observe that to fully describe low energy tttt— scattering 
in the Chiral Resonance framework a broad scalar (t(550) particle is needed. 
In the second part of the thesis we investigate the heavy baryon spectra in the 
bound state picture. In this picture the heavy baryon is treated as a heavy spin 
multiplet of mesons {Qq) bound in the background field of the nucleon (qqq), 
which in turn arises as a soliton configuration of light meson fields. We show 
that a relativistic model with light vectors gives a very satisfactory account of 
the S* — Sc hyperfine splitting in contrast to the model without light vectors. 
In the last chapter we also present a generalization of the bound state model. 
Indeed by binding heavy spin excited multiplets to the background Skyrmion 
field we can describe the spectrum of excited heavy baryons of arbitrary spin. 
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INTRODUCTION 



The non abelian gauge theory which describes, in the perturbative regime, 
the strong interactions is Quantum Chromodynamics (QCD). Quarks and 
gluons are the fundamental degrees of freedom of the theory. A key feature 
of the theory (due to quantum corrections) is asymptotic freedom, i.e. the 
strong couphng constant increases as the energy scale of interest decreases. 
The perturbative approach becomes unreliable below a characteristic scale 
of the theory (A). Quarks and gluons confine themselves into colorless 
particles called hadrons (pions, protons,..). The latter are the true physical 
states of the theory. 

We need to investigate alternative ways to describe strong interactions, 
and in general any asymptotically free theory, in the non perturbative 
regime. This is the fundamental motivation of the present thesis. Al- 
though the underlying gauge theory cannot be easily treated in the non 
perturbative regime we can still use its global symmetries as a guide to 
build Effective Lagrangian Models. These models will be written directly 
in terms of the colorless physical states of the theory, i.e. hadrons. Two 
relevant, approximate symmetries are of extreme theoretical as well as 
phenomenological interest: Chiral Symmetry and Heavy Spin Symmetry. 
Chiral Symmetry is the approximate symmetry associated with the light 
quarks u, d and s, when their masses niq are considered to be negligible 
compared to the invariant scale of the theory A. Heavy Spin Symmetry 
on the other hand is associated with the heavy quarks b, c, and t, whose 
masses Mq are large with respect to A. These symmetries greatly reduce 
the number of unknown parameters and provide a systematic expansion 
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of the effective Lagrangians in rUg and 1/M. Another useful guide for 
the construction of effective Lagrangians will be the expansion for large 
number of colors Nc of QCD. In this thesis we will show that the effective 
Lagrangian models, encoding the relevant symmetries of the underlying 
fundamental gauge theory, provide a reasonable understanding of the in- 
teractions among light mesons at intermediate energies as well as of the 
properties of Heavy baryons. 

In part I of the thesis we will focus our attention on the physics related 
to the three light quark flavors u, d and s. In part II we will consider the 
physics associated with the heavy quarks. 

In the first part of the thesis we will show that it is possible to build a 
reasonable chiral and crossing symmetric model for light meson interac- 
tions at intermediate energies, which we have called the Chiral Resonance 
Model (ChRM). Using the previous scheme we will demonstrate that we 
can understand tttt— scattering up to the 1 GeV region by considering 
the resonance exchange together with the contact term contributions and 
by employing a suitable regularization procedure. We will also observe 
that to fully describe low energy tttt— scattering in the Chiral Resonance 
framework a broad scalar (t(550) particle is needed. At the time when 
the research work was done this state was not present in the Particle data 
Group review (PDG). In the latest PDG this state is finally present. Al- 
though the parameters associated with the a are still not well known the 
quoted range of parameters (mass and width) is consistent with the one 
determined using the Chiral Resonance Model. We will also see that the 
resonance /o(980) has to be introduced with a mechanism a la Ramsauer- 
Townsend. It seems likely that any crossing symmetric approximation will 

XV 



have a similar form. We will regard the ChRM as a leading order 1/Nc 
mean field approximation for Quantum Chromodynamics. 

In the second part of the thesis we will investigate the heavy baryon spec- 
tra in the bound state picture. In this picture the heavy baryon is treated 
as a heavy spin multiplet of mesons {Qq) bound in the background field of 
the nucleon {qqq) , which in turn arises as a soliton configuration of light 
meson fields. A nice feature of this approach is that it permits, in prin- 
ciple, an exact expansion of the heavy baryon properties in simultaneous 
powers of 1/M and l/Nc, where M is the heavy quark mass. We will 
show that a relativistic model with light vectors gives a very satisfactory 
account of the E* — Ec hyperfine splitting in contrast to the model without 
hght vectors. We will also show that the source of hyperfine splitting is 
hidden in non manifest heavy spin breaking terms. In the last chapter we 
will present a generalization of the bound state model. Indeed by binding 
heavy spin excited multiplets to the background Skyrmion field we will see 
that we can describe the spectrum of excited heavy baryons of arbitrary 
spin. 
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Chiral Resonance Model 



1 



Chapter 1 

Introduction to the Chiral 
Resonance Model 

1.1 Brief review of QCD, Chiral Symmetry and 



The non abelian theory which describes, in the pcrturbative regime, the strong inter- 
actions is Quantum Chromdynamics (QCD). Quarks and gluons are the associated 
degrees of freedom [1,2]. SU(3) is the non abehan gauge group of QCD. We will 
denote by the fermionic matter fields associated with quarks, a — u, d, s, ... is the 
flavor index, while c = 1, 2, 3 is the color index. The gauge bosons {gluons) belong to 
the adjoint representation of the gauge group (G^, a = 1, ... ,8). The classical QCD 
Lagrangian density (i.e. which does not include quantum corrections) is 




a 



(1.1) 
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where the color indices are summed via the Einstein convention. The field strength 
tensor is 

^Ixu — 9fiGl — di,G^^ — QsP^'^GbnGcv , (1-2) 

where Qs is the coupling constant and /"^"^ are the group structure constants. The 
covariant derivative acting on a single quark is 

D,q^{d, + igsGlT'^)q, (1.3) 

where T" are the SU(3) generators. A key feature of the theory, due to quantum 
corrections, is asymptotic freedom, i.e. the strong coupling constant gs, decreases as 
the renormalization scale (//) increases [3-5]. Indeed at first order in perturbation 
theory quantum corrections provide 

= Tx H , (1-4) 

(33 - 27V,) In (^) 

where Nf is the number of fiavors whose masses are less than jji and A is a renormal- 
ization invariant scale while, g^^ = Anag. It is clear that the pcrturbative approach 
is not reliable when — > A. Hence Eq. (1.4) itself cannot be trusted in this regime. 
However, a commonly accepted working hypothesis is that the strong coupling con- 
stant keeps increasing at low energies, leading to the phenomenon of quark and gluon 
confinement. In Eq. (1.4) as is a function of A. So the latter can actually be sub- 
stituted for the dimension free quantity in any result deduced using QCD. This 
phenomenon is called dimensional trasmutation. 

We need to investigate alternative ways to describe strong interactions, and in 
general any asymptotically free theory, in the non perturbative regime. This is the 
fundamental motivation of the present thesis. 

It is, however, possible to deduce more information from QCD, by studying it for 
some limiting values of its parameters. The parameters which we can tune are ag 
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(i.e. A), the flavor number, the quark masses and the color number ^. 

For the purposes of the research work presented in part I we will consider only the 
lightest three flavors u, d and s {Nf — S). In part II we will consider heavy quark 
physics. The QCD Lagrangian density in the limit of zero mass for the light quarks 
gains the classical symmetry Ul{3) ® Ur{3). The quark fields will transform as (for 
convenience we omit the color indices) 

ql^gL^L: Q%^94Q^: (1-5) 

where g^^R G ?7l,_r(3) acts, respectively, on the Left and Right components of the 
Dirac spinors. 

The classical symmetry C/l(3)®C/r(3) is explicitly broken via quantum corrections, 

i.e. 

Ul{S) ® Ur{3) Uv{l) ® SUl{3) ® SUr{3) , (1.6) 

This phenomenon is called the chiral anomaly. Uv{l) is associated with baryon num- 
ber conservation. Invariance under the non abelian symmetry group SUl{3)<^SUr{3) 
predicts that all the physical states must be classified according to irreducible repre- 
sentations of the group. Nature docs not realize such a physical spectrum and it is 
assumed that the symmetry is spontaneously broken via 

SUl{3) ® SUr{3) ^ SUv{3) . (1.7) 

The Nambu-Goldstone theorem states that the lost symmetry must be compensated 

by the presence of massless bosons (Goldstone bosons). In QCD such Goldstone 

bosons, associated with the spontaneously broken chiral symmetry, can be identified 
^ Among these parameters (for completness) we must mention the O — vacuum parameter, asso- 
ciated with the strong CP (Charge times Parity) violation. The experimental determination of the 
neutron dipole moment provides a very stringent upper bound, i.e. < 2 x 10~^°. Hence, in the 
following discussion we will not consider it. 
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with the octet of pseudoscalar mesons tt, K, K, rj. All of the other particles can 
be classified according to the irreducible representations of the diagonal subgroup 
SUv{3). 

Now we will consider the large limit of Quantum Chromodynamics [6]. The 
expansion in 1/Nc will allow us to deduce much new information whose key feature 
is its validity in the non perturbative as well as perturbative regime of QCD. In this 
paragraph we will review only some of the 1/Nc predictions, the reader will find more 
results in [6]. In particular here we will consider only the key features relevant for 
part 1 of the thesis. 

For Nc —>■ oo, keeping constant the following product agNc one can prove the 
following meson properties. 

a At the leading order in l/Nc the two point Green's functions built out of local 
quark bilinear operators (J= qq, qjfj,q, etc.) are saturated by the exchange of 
an infinite number of non exotic meson resonances, i.e. of the qq type. 

<Ak)A-k)>=j:—^. (1.8) 

rrin is n}^ meson mass, and a„ =< 0| J|n > is the J matrix element which creates 
the -n}^ meson from the vacuum. 

b The meson-meson scattering amplitude described by a 4-point Green's function, 
schematically represented in Fig.l.l(a), is 0{1/N^. So we can deduce that the 
interactions among the mesons are subleading in the expansion. 

c Non exotic mesons are stable. The 3-point Green function associated with the 
decay amplitude is represented in Fig. 1.1(b) and is 0{l/\/W^. Hence the decay 
itself is 0{l/N^). 

The chiral limit together with the 1/A^c expansion are the two ingredients moti- 
vating the chiral resonance model which we will describe in the following paragraph. 
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a b 

Figure 1.1: (a) schematically represents the 4-point Green function, (b) 
represents the 3-point Green function. In (a) and (b) the external lines 
represent a quark loop, while (•) are the insertions of the non exotic meson 
fields. 

1.2 Chiral Resonance Model 

Asymptotic freedom strongly reduces the predictive power of QCD, in particular for 
the description of strong interactions among mesons at low energies. Chiral perturba- 
tion theory (xPT) [7] allowes a systematic investigation of the strong interactions of 
the octect of Goldstone bosons (tt, K, K, t]). xPT, essentially, relies on the expansion 
in energies of the scattering amplitudes, and it improves the Born terms of the Chiral 
Lagrangian by including quantum corrections via loops and necessary counterterms. 
The energy expansion converges very fast for energies less than 400 — 500 MeV. It 
is a very hard task to extend the xPT scheme at higher energies, since the physical 
singularities associated with the resonances which exist in this region cannot be easily 
reproduced by a truncated (in energy) power expansion. 

In order to describe the interactions among light mesons up to an energy (in the 
center of mass) of 1 — 1.5 GeV is clear that one is forced to include the effects of these 
resonances in this region. 

The question is: How to include such resonances ? As noted in the previous 
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paragraph the meson-meson scattering amphtudes in the leading order in 1/Nc can 
be obtained by summing up all possible tree diagrams obtained via an effective La- 
grangian density, which contains an infinite number of bosonic resonances per each 
possible spin. An infinite number of contact terms is also allowed [6] . 

The 1/Nc analysis is very interesting, since it supports the Lagrangian effective 
models as reasonable models, but at the moment is so general as to appear useless. 
We will see that this is not the case. 

Using the l/Nc analysis as a guide we will try to build a new effective model [8~10] 
to describe the meson scattering up to about 1 GeV. 

A scattering amplitude built using an effective Lagrangian, automatically satisfies 
crossing symmetry. On the other hand just calculating the tree approximation to an 
effective Lagrangian will not guarantee that unitarity is satisfied. This is the handle 
we will use to try to investigate additional structure. Unitarity has of course the 
consequence that the amplitude must have some suitable imaginary term which in 
the usual field theory is provided by loop diagrams. However the leading l/Nc ap- 
proximation will give a purely real amplitude away from the singularities at the direct 
s — channel poles. We may consider the imaginary part of the leading amplitude 
to consist just of the sum of delta functions at each such singularity. Clearly, the real 
part has a much more interesting structure and we will mainly confine our attention 
to it. 

Unitarity has the further consequence that the real parts of the partial wave 
amplitudes must satisfy certain well known bounds. The crucial question is how 
these bounds are satisfied since, as will see, individual contributions tend to violate 
them badly. At first one might expect that all of the infinite number of resonances are 
really needed to obtain cancellations. However the success of chiral dynamics at very 
low energies where none of the resonances have been taken into account suggests that 
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this might not be the case. It is clear that we need to postulate a new principle which 
has been introduced in Ref. [8] and we will call it local cancellation. Local cancellation 
cannot be easily derived from QCD, but if true it would greatly simplify the task of 
extending the phenomenological description of scattering processes to higher energies. 

According to this simple principle the inclusion of the resonances, up to the en- 
ergy one wants to describe, in the direct channel as well as in the cross channels, 
will allow us to saturate the unitarity bounds and also will allow us to describe the 
phenomenology involved. Now one can easily understand the relevant role played by 
crossing symmetry, already pointed out in Veneziano's famous paper [11]. Indeed, as 
we will see in the next chapter, local cancellations among the different contributions 
from the resonances will, non trivially, enforce unitarity. In tttt scattering, we will 
see that the introduction of the cross term due to p meson exchange in a chirally 
invariant manner, in the Isospin = and zero orbital angular momentum case, sub- 
stantially delays the onset of the severe unitarity violation which would be present in 
the simplest chiral Lagrangian of pions. 

The onset of the unitarity bound via the local cancellation principle will also 
require the existence of the scalar particle (t(550) [8-10] in order to succesfuUy describe 
TTTT— scattering up to about 1.2 GeV. 

In order to start this investigation it seems reasonable to keep in the Lagrangian 
only local operators with the lowest number of derivatives. Operators with a higher 
number of derivatives destabilize the theory at high energy and will make it more 
difficult to accomplish the required local cancellation. Further restrictions will be 
obtained by imposing chiral symmetry. 

In the last paragraph of the present chapter we will show how to include the scalar, 
vectorial and tensorial resonances in order to preserve chiral symmetry. 

We will now summarize the chiral resonance model according to the following 
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logical scheme. 

a The effective Lagrangian will contain only terms with the lowest number of 
derivatives. These local operators will be chirally symmetric. 

b The model will be suitable for the real part of the physical amplitude. Such an 
amplitude must obey crossing symmetry. 

c We will include in the scheme only the resonances whose masses lie in the energy 
range of interest. Away from the poles in the direct channel unitarity is expected 
to be provided by the local cancellation principle. 

d A suitable regularization method must be employed to avoid divergences at each 
direct channel pole for the physical amplitude. 

In the next paragraph we will explain how to regularize the s— channel poles. 
The other points will be investigated in the next two chapters. In particular we will 
study the tttt —>■ tttt and tttt KK scattering amplitudes. These two channels are a 
classically important check for any model whose goal is to describe meson interactions 
at intermediate energies. 



1.3 Regularization of the Model 

In the large Nc picture the leading amplitude (of order is a sum of polynomial 

contact terms and tree-type resonance exchanges. Furthermore the resonances should 
be of the simple qq type; glueball and multi-quark meson resonances are suppressed. 
In our phenomcnological model there is no way of knowing a priori whether a given 
experimental state is actually of the qq type. For definiteness we will keep all relevant 
resonances even though the status of a low lying scalar resonances hke the /o(980) 
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has been considered especially controversial [32]. If such resonances turn out in the 
future to be not of type qq, their tree contributions would be of higher order than 
l/Nc- In this event the amplitude would still of course satisfy crossing symmetry. 

The most problematic feature involved in comparing the leading l/Nc amplitude 
with experiment is that it does not satisfy unitarity. In fact, resonance poles like 

will yield a purely real amplitude, except at the singularity, where they will diverge 
and drastically violate the unitarity bound. Thus in order to compare the l/Nc 
amplitude with experiment we must regularize the denominators in some way. The 
usual method, as employed in Ref. [8], is to regularize the propagator so that the 
resulting partial wave amplitude has the locally unitary form 

. (1.10) 

M^-s- iMT ^ ' 

This is only valid for a narrow resonance in a region where the background is negligible. 
Note that the —iMT is strictly speaking a higher order in l/Nc effect. 

For a very broad resonance there is no guarantee that such a form is correct. 
Actually, in Ref. [8] it was found necessary to include a rather broad low lying scalar 
resonance (denoted a (550)) to avoid violating the unitarity bound. A suitable form 
turned out to be of the type 

MG 

(1.11) 



iVP -s- iMG' ' 

where G is not equal to the parameter G' which was introduced to regularize the 
propagator. Here G is the quantity related to the squared coupling constant. 

Even if the resonance is narrow, the effect of the background may be rather im- 
portant. Demanding local unitarity in this case yields a partial wave amplitude of 
the well known form [13]: 



e 



2ii 



''Mr 



+ e^^sin(5, (1.12) 



M2 - s - iMV 
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where 5 is a background phase (assumed to be slowly varying) . We will see in the next 
chapters that such a regularization method is needed to fully understand the tttt — > tttt 
and TTTT — > KK scattering in the /o(980) region. We will adopt a point of view in 
which this form is regarded as a kind of regularization of our model. Of course, non 
zero d represents a rcscattcring effect which is of higher order in l/Nc- The quantity 
e^*'^, taking 5 — constant, can be incorporated, for example, into the squared coupling 
constant connecting the resonance to two pions. In this way, crossing symmetry can 
be preserved. Prom its origin, it is clear that the complex residue does not signify the 
existence of a ghost particle. The non-pole background term in Eq. (1.12) and hence 
S is to be predicted by the other pieces in the effective Lagrangian. 

Another point which must be addressed in comparing the leading 1/Nc amplitude 
with experiment is that it is purely real away from the singularities. The regulariza- 
tions mentioned above do introduce some imaginary pieces but these are clearly very 
model dependent. Thus it seems reasonable to compare the real part of our predicted 
amplitude with the real part of the experimental amplitude. Note that the difficulties 
mentioned above arise only for the direct channel poles; the crossed channel poles 
and contact terms will give purely real finite contributions. 

It should be noted that if we predict the real part of the amplitude, the imaginary 
part can always be recovered by assuming elastic unitarity. 

1.4 Chiral Lagrangian 

In the low energy physics of hadrons it is important to correctly introduce the spon- 
taneous chiral symmetry breaking structure. In our approach we will introduce it 
via non linear realizations. It is known that this method reproduces the low energy 
theorems obtained via current algebra. We start here with the definition of the 3 x 
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3 matrix U, 

, e = e"^/^^ , (1.13) 

where is a pion decay constant. U is parameterized by the pseudoscalar matrix 
4>, which is identified with the pseudoscalar meson octect. Under the chiral group 
C/(3)l X C/(3)r, U transforms hnearly. 

U ^ 9LUgi , (1.14) 

where 5fL,R ^ U(3)lj^. Under the chiral transformation Eq. (1.14), ^ transforms non- 
linear ly: 

e ^ g^CK^{<f>,9L,9R) = K{<f>,g^,g^)C9l . (1.15) 

The previous equation implicitly defines the matrix K{(f), gj^, g^i). The vector meson 
nonet is introduced as a gauge field [15] which transforms as 

^ Xp^Xt _ lx9^xt ^ (1.16) 

where ^ is a gauge coupling constant. (For an alternative approach see, for a review, 
Ref. [16].) It is convenient to define 



{^d,e + ^^d,^) , (1.17) 



which transform as 



V 



Kv^K'' +iKdi,K'' , (1.18) 



and are the chiral group Maurer — Cartan one form. Using the above quantities we 
construct the chiral Lagrangian including both pseudoscalar and vector mesons: 
1 



{gp, + v,Y\ + -fTv \p^p^] - -lY [F,,{p)F^^{p)] , (1.19) 
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where F^iy — d^j,pv — dyp^ + ig\p,x., pv\ is a gauge field strength of vector mesons, g is 



connected to the vector meson couphng to two pions gp^^T^ via 

ml 

gpnir 



[1.20) 



In the real world chiral symmetry is explicitly broken by the quark mass term —fnqAiq, 
where m = (m^ + md)/2, and is the dimensionless matrix: 



M 



I 



\ 



1 + y 



X 



:i.2i) 



Here x and y are the quark mass ratios: 



X 



m 



y 



1 / mg - rriu 

2 V m 



(1.22) 



These quark masses lead to mass terms for pseudoscalar mesons. Moreover, in con- 
sidering the processes related to the kaon, (in this thesis we will consider tttt — > KK 
scattering amplitude) we need to take account of the large mass gap between the s 
quark mass and the u and d quark masses. These effects are included as SU(3) sym- 
metry breaking terms in the above Lagrangian, which are summarized, for example, 
in Refs. [17, 18]. Here we write the lowest order pseudoscalar mass term only: 



f A'Tr 

mass u ±i 



:i.23) 



where 5' is an arbitrary constant. 

We next introduce higher resonances into our Lagrangian. Firstly, we write the 
interaction between the scalar nonet field S and pseudoscalar mesons. Under the 
chiral transformation, this S transforms as — > KSK^ . A possible form which 
includes the minimum number of derivatives is proportional to 



(1.24) 
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The coupling of a physical isosinglet field to two pions is then described by 




7o 



(1.25) 



Here we should note that chiral symmetry requires derivative-type interactions be- 
tween scalar fields and pseudoscalar mesons. Secondly, we represent the tensor nonet 
field by T^i, (satisfying T^j, — Tj^^, and Tj^ — 0), which transforms as 



The heavier vector resonances such as p(1450) can be introduced in the same way as 
p in Eq. (1.19). 

In general, for a complete description of the mesonic processes one should include 
the effects due to chiral anomahes. Since in this part of the thesis we will not deal 
with such processes we refer the reader to the following references [17-19]. 



T, 




(1.26) 



The interaction term is given by 



(1.27) 



Chapter 2 



A fundamental process: tttt 
scattering 

2.1 Introduction to tttt scattering 

Historically, the analysis of tttt scattering has been considered an important test of 
our understanding of strong interaction physics (QCD, now) at low energies. It is 
commonly accepted that the key feature is the approximate spontaneous breaking of 
chiral symmetry. Of course, the kinematical requirements of unitarity and crossing 
symmetry should be respected. The chiral perturbation scheme (xPT) [7] , which im- 
proves the tree Lagrangian approach by including loop corrections and counterterms, 
can provide a description of the scattering up to the energy region slightly above 
threshold (400 — 500 MeV). Hence xPT is very useful for describing kaon decays into 
27r and Stt (for an apphcation see Ref. [20]). In Ref. [21] the derivative expansion has 
been used to estimate some of the parameters associated with CP violation for the Stt 
decays of neutral kaons. In order to describe the scattering up to energies beyond this 
region (say to around 1 GeV) it is clear that the effects of particles lying in this region 
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must be included. The Chiral Resonance Model (ChRM) provides a simple prescrip- 
tion for including these resonances. We will see that by applying the ChRM recipe 
we will be able to understand the scattering processes schematically represented in 
Fig. 2.1. The kinematics, the unregularized invariant scattering amplitudes and the 



Figure 2.1: Schematic representation of a scattering process. 



generic scattering matrix S are presented in Appendix A. 

In the previous chapter we stated that the ChRM predicts the real part of the 
scattering amplitude. The imaginary term can be recovered via the unitarity relations. 
Specializing Eq. (A. 6) in Appendix A to the tttt channel we have, for the imaginary 
piece // of the /, / partial wave amplitude 



^'4 



1 ± ^n{' - AR{^ 



(2.1) 



where tjI is the elasticity parameter. / and / are the isospin and the orbital angular 
momentum. Obviously, this formula is only meaningful if the real part obeys the 
bound 

\Rl\ < I . (2.2) 

The main difficulty one has to overcome in obtaining a unitary amplitude by the 
present method is the satisfaction of this bound Eq. (2.2). 

We want to remark that making the rcgularizations of direct channel poles such as 
in Eqs. (1.10) and (1.12) which provide unitarity in the immediate region of a narrow 
resonance, is not at all tantamount to unitarizing the model by hand. 
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In this chapter we will present the first evidence of local cancellation, obtained 
among the current algebra and the vector meson p contributions. We will then show 
that we need to introduce a broad scalar resonance (i.e. F > M) to fully understand 
the scattering process at low energies (< 800 MeV). The effect of chirally invariant 
contact terms with a higher number of derivatives will also be investigated. 

2.2 Current Algebra and p vector— meson exchange 

In this section we will study the partial waves for tttt scattering computed in a chiral 
Lagrangian model which contains both the pseudoscalar and vector mesons, (i.e., the 
lowest lying s-wave quark-antiquark bound states). Near threshold tjI — 1, Rj is 
small and we should choose the minus sign in Eq. (2.1) so that 



In the large Nc limit the amplitude near threshold is purely real and of the order 
l/Nc- This is consistent with Eq. (2.3) which shows that //(s) is of order l/N^ and 
hence comes in at the second order. This agrees with the chiral perturbation theory 
approach [7] in which Rj [s) comes from the lowest order tree diagram while // arises 
from the next order loop diagram. On the other hand, when we depart from the 
threshold region the 1/Nc approach treats the contribution of the p-meson at first 
order while the chiral perturbation theory approach treats it at second and higher 
orders. A straightforward computation using the pion lagrangian (second term in 
Eq. (1.19)) together with the explicit chiral breaking term in Eq. (1.23) yields the tttt 
scattering amplitude [22] defined in Eq. (A. 8): 



lUs) ^ [Rif. 



(2.3) 



AcA{s,t,u)^2 



s — 




(2.4) 
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This equation will be called the current algebra result. With (A. 9) and (A. 11) we 
obtain Rq{s) — T'{\.q(s) as illustrated in Fig. 2.2. The experimental Roy curves [23] 
are also shown. Up until about 0.5 GeV the agreement is quite reasonable (and can 
be fine-tuned with second order chiral perturbation terms) but beyond this point Rq 
keeps increasing monotonically and badly violates the unitarity bound (2.2). We 

1 .00 
0.75 
0.50 
0.25 

0.00 

0.25 0.35 0.45 0.55 0.65 0.75 

Vs GeV 

Figure 2.2: The solid line is the current algebra result for Rq. The dotted 
and dot-dashed lines are the Roy curves for Rq. 

will see now that the introduction of the p-meson greatly improves the situation. The 
chiral Lagrangian for the vector-meson pion system is displayed in Eq. (1.19). The 
Lagrangian piece in (1.19) yields both a pole- type contribution (from the P/^v^ cross 
term) and a contact term contribution (from the v^v^ term) to the amplitude at tree 
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level: 



Ms:t,u) = -^-^^{^^ + -^]+^f^[{t-s) + {u-s)] . (2.5) 
2 \mj — t mj — uj 2mj 

where rUp — 0.769 GeV and g'p^^ = 8.56. We notice that the entire first term in (1.19) 



Tr 



{~gp^. + v^f , (2.6) 



is chiral invariant since V/^ and —gp^ transform identically. However the Tr(p^^;'*) and 
Tr(f^f^) pieces are not separately chiral invariant. This shows that the addition of 
the p-meson in a chiral invariant manner necessarily introduces a contact term in 
addition to the minimal pole term. Adding up the two terms in Eq. (2.5) and the 
term in Eq. (2.4) yields finally 



.2 ^2 

AcA+p{s,t,u) = 2 



t{u — s) u(t — s) 
— t ml — u 



(2.7) 



F2 2ml 

In this form we see that the threshold (current algebra) results are unaffected since 
the second term drops out at t = m = 0. An alternative approach [24] to obtaining 
Eq. (2.7) involves introducing a chiral invariant pTrvr interaction with two more deriva- 
tives. ^4(5, t, u) has no singularities in the physical region. Reference to Eq. (A. 9) 
shows that the isospin amphtudes T{\ and Tp^ also have no singularities. However 
the T^-^ amplitude has the expected singularity at s = m^. This may be cured in a 
conventional way, while still maintaining crossing symmetry, by the replacements 

1 1 



m^ — t,u m^ — t,u — impTp6{t, u — i.m^) 
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A modification of this sort would enter automatically if we were to carry the computa- 
tion to order — -. However we shall regard (2.8) as a phenomenological regularization 
of the leading amplitude. 

Now let us look at the actual behaviour of the real parts of the partial wave 
amplitudes. i?Q, as obtained from Eq. (2.7), is graphed in Fig. 2.3 for an extensive 
range of -y^, together with the current algebra result. 
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We immediately see that there is a remarkable improvement; the effect of adding 
p is to bend back the rising Rq{s) so there is no longer a drastic violation of the 
unitarity bound until after ^/s — 2 GeV. There is still a relatively small violation 
which we will discuss later. Note that the modification in Eq. (2.8) plays no role in 
the improvement since it is only the non-singular t and u channel exchange diagrams 
which contribute. 

It is easy to see that the c?e/a?/erf drastic violation of the unitarity bound \Rj\ < - is 
a property of all partial waves. We have already learned from (2.7) that the amphtude 
A{^s^t^u) starts out rising linearly with s. Now Eq. (2.5) shows (for fixed scattering 
angle) that for large s the p exchange terms behave as The leading large s behavior 
will therefore come from the sum of the original current- algebra term and the new 
contact-term: 



But k is numerically around 2 [25], so A{s,t,u) eventually decreases linearly with s. 
This turn-around, which is due to the contact term that enforces chiral symmetry, 
delays the onset of drastic unitarity violation until well after the p mass. It thus seems 
natural to speculate that, as we go up in energy, the leading tree contributions from 
the resonances we encounter (including both crossed channel as well as s-channel 
exchange) conspire to keep the -Rf(s) within the unitarity bound This is the first 
evidence of the local cancellation "phenomenon. We notice that the cancellation among 
resonance and contact term contributions for large s allowed Weinberg in Ref. [26] to 
deduce some remarkable asymptotic relations about the resonance spectrum and the 
pion couplings. 

In Figure 2.4 we show the partial waves R\ and ll computed using Eq. (2.7) 
and Eq. (2.8). Not surprisingly, these display the standard resonant forms. The 
dominance of the vector-meson p in the I — I — 1 channel is used in xPT to estimate 




(2.9) 
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some of the counterterm coefficients [27]. For completeness we present the R^ and R2 
amphtudes in Fig. 2.5. We notice that these amphtudes obey the unitarity hmits up 
to about 2 GeV. 

In this section we have shown that the inclusion of the p-meson dramatically re- 
duces the onset of unitarity violation in the Rq channel (see. Fig. 2.3). We notice that 
this channel contains the exchange vector meson diagram which provides a decisive 
contribution. This result is a confirmation, not only, of the exchange symmetry but 
also of the relevance of the crossing symmetry. 

The cancellation for large s is due to the contact term v^Vfj,. The last term is 
fundamental in order to introduce the p in a chiral symmetric way. In the following 
we will mainly study the Rq channel, since this has the worst behaviour for large s. 

2.3 Is the cr(550) alive ? 

As shown in Figure 2.6, although the introduction of the p-meson dramatically im- 
proves unitarity up to about 2 GeV, Rq violates unitarity to a lesser extent starting 
at around 500 MeV. What is needed to restore unitarity over the full range of interest 
and to give better agreement with the experimental data for ^/s < 900 MeVl 

i. Below 450 MeV, Rq{s) actually hes a little below the Roy curves (see Fig. 2.2). 
Hence it would be nice to find a tree level mechanism which yields a small 
positive addition in this region. 

ii. In the 600 — 1300 MeV range, an increasingly negative contribution is clearly 
required to keep Rq within the unitarity bound. 

In Reference [8] it has been shown that it is possible to satisfy both of these criteria by 
introducing a broad scalar resonance (like the old a) with a mass around 550 MeV. 
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Figure 2.5: The solid line is the n + p contribution for the real part of 
1 = 2, 1 = 0. The dot-dashed line is the tt + p contribution for the real 
part of / = 0, 1 = 2. 



CHAPTER 2. A FUNDAMENTAL PROCESS: tttt SCATTERING 



2.0 




Figure 2.6: Enlarged version of Fig. 2.3. The solid line which shows the 
current algebra + p result for Rq is much closer to the unitarity bound of 
0.5 than the dashed hne which shows the current algebra result alone. 
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In this model the a is not reahzed as a chiral partner hut is introduced as a matter field 
with respect to chiral transformations. The a contribution to the invariant amplitude 
A{s, t, u) is 

where ^ 

and G is related to the couphng constant 70 defined in Eq. (1.25) as 

Note that the factor (s — 2m^)^ is due to the derivative-type coupling required 
for chiral symmetry in Eq. (1.25). The total amplitude will be crossing symmetric 
since ^4(5, t, u) and A{u, t, s) in Eq. (A. 8) are obtained by performing the indicated 
permutations. G' is a parameter which we introduce to regularize the propagator. It 
can be called a width, but it turns out to be rather large so that, after the p and tt 
contributions are taken into account, the partial wave amplitude does not clearly 
display the characteristic resonant behavior. In the most general situation one might 
imagine that G could become complex as in Eq. (1.12) due to higher order in l/N^ 
corrections. It should be noted, however, that Eq. (1.12) expresses nothing more than 
the assumption of unitarity for a narrow resonance and hence should not really be 
applied to the present broad case. 

A reasonable fit was found in Ref. [8] for G purely real, but not equal to G'. By the 
use of Eq. (2.1), unitarity is in fact locally satisfied. In [10] a best overall fit is obtained 
with the parameter choices M„ = 559 MeV, G/G' = 0.29 and G' = 370 MeV. These 
have been slightly fine-tuned from the values in Ref. [8] in order to obtain a better 
fit in the 1 GeV region. The contribution to due only to the presence of the low 
mass broad scalar (t(550) is displayed in Fig. 2.7. 
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Figure 2.7: Contribution to Rq due to the (j(550) particle. 
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The curve crosses zero at about ^/s = 560 MeV. The small deviation is due to the 
effect of the cross diagrams. The result for the real part of Rq due to the inclusion 
of the a contribution along with the tt and p contributions is shown in Fig. 2.8. It 
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Figure 2.8: The sohd line is the current algebra + p + a result for Rq. 
The experimental points, in this and suceeding figures, are extracted from 
the phase shifts using Eq. (A. 6) and actually correspond to i?Q/?7Q. (□) 
are extracted from the data of Ref. [28] while (A) are extracted from the 
data of Ref. [29]. The predicted Rq is small around the 1 GeV region. 

is seen that the unitarity bound is satisfied and there is a reasonable agreement with 
the experimental points [28, 29] up to about 800 MeV. Beyond this point the effects 
of other resonances (mainly the /o(980)) arc required. From Eqs. (2.10), (A. 9) and 
(A. 11) we see that the contribution of a to Rq becomes negative when s > M^. This 
is the mechanism which leads to satisfaction of the unitarity bound (c.f. Fig. 2.6). 
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For s < one gets a positive contribution to Rq. This is helpful to push the 
predicted curve upwards and closer to the experimental results in this region, as 
shown in Fig. 2.9. The four-derivative contribution in the chiral perturbation theory 

0.60 
0.50 
0.40 
0.30 
0.20 
0.10 
0.00 

0.20 0.30 0.40 0.50 0.60 

\|T (GeV) 

Figure 2.9: A blowup of the low energy region. The sohd line is the 
current algebra + p contribution to Rq. The dashed line includes the a 
particle and has the effect of turning the curve down to avoid unitarity 
violation while boosting it at lower energies. 

approach performs the same function, however it does not change sign and hence does 
not satisfy the unitarity bound above the 450 MeV region [30] . 

It is also interesting to notice that the main effect of the sigma particle comes 
from its tail in Fig. 2.7. Near the pole region, its effect is hidden by the dominant 
TT + p contribution. This provides a possible explanation of why such a state may 
have escaped definitive identification. It is interesting to remark that a particle with 
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mass and width very similar to those given above for the a was predicted [31] as 
part of a multiquark qqqq nonet on the basis of the MIT bag model. Hence, even 
though they do not give rise to formally leading tttt amplitudes in the 1/Nc scheme, 
the picture has a good deal of plausibility from a polology point of view. It is not 
hard to imagine that some l/N,. sublcading effects might be important at low energies 
where the QCD coupling constant is strongest. For completeness in Fig. 2.10 we also 
show the real part of the 1 — 21 — scattering amphtude (Rl). We have used 
for the a parameters those presented in the first column of Table 3.2. It is clear 
that there is a fair agreement with the data up until about 1 GeV for the dashed line 
which represents the current algebra +p+ a, while the current algebra prediction alone 
departs quite soon from the experimental results. The solid line {current algebra +p) 
seems to better follow the experimental data up to about 1.2 GeV. However since 
we have not tried to fit this channel we will merely stress that the results are also 
consistent with the local cancellation principle. For other authors [26, 32, 35] the a is 
a qq state, while in [36] is considered as an interpolating multi pionic state. In the 
actual model it is not possible to uncover what makes the a. This state is however 
essential according to the Chiral Resonance Model. M^., G and G' are in practice 
the only unknown parameters in the present model. Here we remark that the a has 
been introduced to understand low energy tttt scattering in the following two papers: 
Exploring tttt Scattering in the 1/Nc Picture (Ref. [8]) and Simple description of tttt 
scattering to 1 GeV {Rei. [10]). At the time when these papers were written this state 
was not present in the Particle Data Group (PDG) [37] review. In the latest PDG [38] 
this state is finally present. Although the parameters associated with the a arc still 
not well known the quoted range of parameters (mass and width) is consistent with 
the one determined in the Chiral Resonance framework. As the conclusion to this 
paragraph we can say that the cr is alive and well. 
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Figure 2.10: The solid line is the current algebra +p contribution for R^. 
The dashed line is the current algebra +p + cr contribution for R^. The 
dot-dashed line is the current algebra contribution for R^. 
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2.4 Comment on higher derivative contact terms 

In the previous paragraph we concluded that the a is crucial in order to satisfy the 
unitarity bounds and to fit the low energy data in the Rq channel. But why do we need 
to introduce a new particle if we still have at our disposal higher derivative contact 
terms ? Here we will analyze this possibility. For simplicity, here we will investigate 
four-derivative contact terms. There are two four-derivative chiral invariant contact 
interactions which are single traces in flavor space: 



oTr 



d^UdM^d^'Ud'^U^] + 6Tr \di,Ud''U^d^Ud''U^] , (2.13) 



where a and b are real constants. The single traces should be leading in the 1/Nc 
expansion. Notice that the magnitudes of a and b will differ from those in the chiral 
perturbation theory approach [7] since the latter essentially also includes the effects 
of expanding the p exchange amplitude up to order s^. The four pion terms which 
result from Eq. (2.13) are: 

L4 = [2a (d^TT ■ d^Trf + {b-a) {d^TT ■ d'^iff] +■■■ . (2.14) 

TT 

This leads to the following contribution to the tttt amplitude: 
1 f\ 

A,{s,t,u) ^-^{a [{t - 2mlf + {u - 2mlf] + {b - a){s - 2mlf} . (2.15) 

TT 

It is reasonable to require that Eq. (2.15) yields no correction at threshold, i.e. at 
s — 4:m'^, t — u — 0. This gives the condition b — —a and leaves the single parameter a 
to play with. In Figure 2.11 we show Rq, as gotten by adding the piece obtained from 
Eq. (2.15) for several values of a to the contribution of 7r-|-p. For a = -1-1.0 x 10~^ the 
four-derivative contact term can pull the curve for Rq down to avoid violation of the 
unitarity bound until around ^/s — 1.0 GeV . The price to be paid is that Rq decreases 
very rapidly beyond this point. We consider this to be an undesirable feature since it 
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would make a possible local cancellation scheme very unstable. Another drawback of 
the four-derivative contact term scheme is that it lowers Ro{s) just above threshold, 
taking it further away from the Roy curves. Hence in the following we will stick with 
the deus ex machina a particle and will keep only the minimal number of derivatives. 




\1T (GeV) 

Figure 2.11: Four-derivative contact term contribution for R^. The solid 
hne corresponds to a = -1-1.0. The dotted hne corresponds to a = -1-0.7. 
The dot-dashed line corresponds to a = 0.5. a is represented in units of 
10-3. 



Chapter 3 



Exploring the 1 GeV region 

3.1 What is happening in the 1 GeV region 

In the previous chapter we have shown [8] that the Chiral Resonance Model is capable 
of describing nn scattering up to about 800 MeV in the I — I — channel (see Fig. 2. 8). 
In this chapter wc will investigate the 800 — 1200 MeV energy range [10]. We will 
study in some detail the /o(980) properties and the problems connected with inelastic 
effects induced by the opening of the KK threshold. 

The neutral resonances which can contribute to the I — I — channel have 
the quantum numbers J^^ = 0+''", 1 and 2++. In the quark model these states 
are naturally interpreted as s and p wave qq bound states. We show in Table 3.1 
the specific ones which are included, together with their masses and widths, where 
available from the Particle Data Group (PDG) [37] hstings. 

Reference to Fig. 2.8 shows that the experimental data for lie considerably 
below the n + p + a contribution between 0.9 and 1.0 GeV and then quickly reverse 
sign above this point. We will now see that this distinctive shape is almost completely 
explained by the inclusion of the relatively narrow scalar resonance /o(980) in a 
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M{MeV) 




Br (27r)% 


(7(550) 


0+(0++) 


559 


370 




Pirn 


1+(1-) 


769.9 


151.2 


100 


/o(980) 


0+(0++) 


980 


40-400 


78.1 


/2(1270) 


0+(2++) 


1275 


185 


84.9 


/o(1300) 


0+(0++) 


1000-1500 


150-400 


93.6 


p(1450) 


1+(1-) 


1465 


310 


seen 



Table 3.1: Resonances included in the tttt — > tttt channel as listed in the 
PDG. Note that the a was not present in the 1994 PDG and is not being 
described exactly as a Breit-Wigner shape; we listed the fitted parameters 
shown in column 1 of Table 3.2 where G' is the analog of the Breit- Wigner 
width. 

suitable manner. One can understand what is going on very simply by starting from 
the real part of Eq. (1.12): 



,_(M2 -s)cos(25) -Mrsin(25) 1 , 



(3.1) 



This expresses nothing more than the restriction of local unitarity in the case of 
a narrow resonance in the presence of a background. We have seen that the diffi- 
culty of comparing the tree level 1/Nc amplitude to experiment is enhanced in the 
neighborhood of a direct channel pole. Hence it is probably most rehable to identify 
the background term -sin(25) with our prediction for R^. In the region of interest, 
Fig. 2.8 shows that Rq is very small so that one expects 5 to be roughly 90° (assuming 
a monotonically increasing phase shift). Hence the first pole term is approximately 

(M^ - s)MV 

(M2 - Sf + M2r2 ' ^ ' 

which contains a crucial reversal of sign compared to the real part of Eq. (1.10). 
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Thus, just below the resonance there is a sudden negative contribution which jumps 
to a positive one above the resonance. This is clearly exactly what is needed to 
bring experiment and theory into agreement up until about 1.2 GeV, as is shown in 
Fig. 3.1. The actual amplitude used for this calculation properly contains the effects 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Ys (GeV) Ys (GeV) 

(a) (b) 

Figure 3.1: (a): The solid hne is the current algebra + p + a + /o(980) 

result for Rq obtained by assuming column 1 in Table 3.2 for the a and 

/o(980) parameters (5r(/o(980) ^ 2n) = 100%). (b): The sohd line is 

the current algebra + p + a + /o(980) result for R^ obtained by assuming 

column 2 in Table 3.2 (Sr(/o(980) ^ 27r) = 78.1%) . 

of the pions' derivative coupling to the /o(980) as in Eq. (2.10). 

The above mechanism, which leads to a sharp dip in the I — J — partial wave 
contribution to the Trvr-scattering cross section, can be identified with the very old 
Rams auer- Towns end effect [40] which concerned the scattering of 0.7 eV electrons 
on rare gas atoms. The dip occurs because the background phase of n/2 causes the 
phase shift to go through tt (rather than 7r/2) at the resonance position. (Of course, 
the cross section is proportional to J2i,j{'^J + 1) sin^((5/).) This simple mechanism 
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seems to be all that is required to understand the main feature of tttt scattering in 
the 1 GeV region. 



3.1.1 The Rams auer- Towns end 'Effect 

Here we will compare the real part of the I — J = partial wave amplitude which 
results from our crossing symmetric model with experimental data. Firstly we will 
consider the sum of the contributions of the current algebra, p-meson, a and /o(980) 
pieces. Then we will add pieces corresponding to the next group of resonances, namely 
the /2(1270), the p(1450) and the /o(1300). In this section we will continue to neglect 
the KK channel. 

The current algebra plus p contribution to the quantity yl(s, t,-u) is defined in 
Eq. (2.7). Note that for the I — J — channel the expression in Eq. (2.7) will yield 
a purely real contribution to the partial wave amplitude. The contribution of the 
low-lying a meson was given in Eq. (2.10). For the important /o(980) piece we have 



ReAfQ(^QS0){s,t,u) = Re 



(3.3) 



_m\ - s- imf^FtotifoWis - ^ml) 

where 5 is a background phase parameter and the real coupling constant 7/o7r7r is 
related to the /o(980) — >• tttt width by 



r(/o(980) ^ tttt) - ^ ^'"^^ 



327r ruf, \ 7X1% ^ ' 



We will not consider 5 to be a new parameter but shall predict it as being 

1 



2 



sin(25) = Rl{s = m%) , (3.5) 



where i?Q is computed as the sum of the current algebra, p, and sigma pieces. Since 
the KK channel is being neglected, one might want to set the regularization parameter 
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^totifo) in the denominator to r(/o(980) vrvr). We shall try both this possibility as 

11 1 -1 r(/o(980) tttt) 
well as the experimental one — -— ^ 78.l7o. 

tot [Jo) 

A best fit of our parameters to the experimental data results in the curves shown 
in Fig. 3.1 for both choices of branching ratio. Only the three parameters G/G', G' 
and Mo- are essentially free. The others are restricted by experiment. Unfortunately 
the total width rtot(/o) has a large uncertainty; it is claimed by the PDG to lie in 
the 40 — 400 MeV range. Hence this is effectively a new parameter. In addition we 
have considered the precise value of mf^^ to be a parameter for fitting purposes. The 
parameter values for each fit are given in Table 3.2 together with the values. It 
is clear that the fits are good and that the parameters are stable against variation 
of the branching ratio. The predicted background phase is seen to be close to 90° 
in both cases. Note that the fitted width of the /o(980) is near the low end of the 
experimental range. The low-lying sigma has a mass of around 560 AleV and a 
width of about 370 MeV . As explained in section 3, we are not using an exactly 
conventional Breit- Wigner type form for this very broad resonance. The numbers 
characterizing it do however seem reasonably consistent with other determinations 
thereof [32,41,42]. 



3.2 Next group of Resonances 

The local cancellation principle together with crossing works up to 1 GeV in energy. 
What about above 1 GeV ? Does the inclusion of the resonances present in this range 
of energy modifies the results obtained at low energy ? 

In the 1/Nc approximation the qq mesons belong to ideally mixed nonets [6]. 
Following the 1/Nc prescriptions, at the first order, only excited p— wave and radially 
excited states in the s— wave must be included. 
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With Next Group 


No p(1450) 


BR{fo{980) 27r)% 


100 


78.1 


78.1 


78.1 


100 


78.1 


78.1 


78.1 


100 




1 


1 


0.8 


0.6 


1 


1 


0.8 


0.6 


1 


M/o(98o) {MeV) 


987 


989 


990 


993 


991 


992 


993 


998 


992 


Ttot (MeV) 


64.6 


77.1 


75.9 


76.8 


66.7 


77.2 


78.0 


84.0 


64.6 


(MeV) 


559 


557 


557 


556 


537 


537 


535 


533 


525 


G' {MeV) 


370 


371 


380 


395 


422 


412 


426 


451 


467 


G/G' 


0.290 


0.294 


0.294 


0.294 


0.270 


0.277 


0.275 


0.270 


0.263 


6 (deg.) 


85.2 


86.4 


87.6 


89.6 


89.2 


89.7 


91.3 


94.4 


90.4 




2.0 


2.8 


2.7 


3.1 


2.4 


3.2 


3.2 


3.4 


2.5 



Table 3.2: Fitted parameters for different cases of interest. 



The neutral members of the p— wave qq nonets have the quantum numbers J^*^ = 
O"*""*", l"*""*", l+~ and 2++. Of course, the neutral members of the radially excited 
s — wave qq nonets have J^'-' = and 1 . Only members of the 0"''+, 1 and 
2++ nonets can couple to two pseudoscalars ^. By G-parity conservation we finally 
note that it is the 7 = member of the O^"*" and 2++ nonets and the I — 1 member of 
the 1 nonet which can couple to two pions. Are there good experimental candidates 
for these three particles? 

The cleanest case is the lighter / = member of the 2~^~^ nonet: the /2(1270) has, 
according to the August 1994 Review of Particle Properties (PDG) [37], the right 
quantum numbers, a mass of 1275 ± 5 MeV, a width of 185 ± 20 MeV, a branching 
ratio of 85% into two pions, and a branching ratio of only 5% into KK. On the other 
hand the /2(1525) has a 1% branching ratio into tttt and a 71% branching ratio into 

^It is possible to write down a two point mixing interaction between 0~+ and radially excited 
particles etc., but we shall neglect such effects here. 
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KK. It seems reasonable to approximate the 2++ nonet as an ideally mixed one and 
to regard the /2(1270) as its non-strange member. 

The p(1450) is the lightest hsted [37] particle which is a candidate for a radial 
excitation of the usual p(770). It has a less than 1% branching ratio into KK but the 
TTTT branching ratio, while presumably dominant, is not yet known. In the following 
analysis, for definitness, we will assume that the p(1450) mainly decays into tttt. 
However, we notice that the ir*(1410), which presumably belongs to the same p(1450) 
SU{3) multiplet has a Kn branching ratio of only 7%. Hence we might also expect a 
small coupling to the two pions. We will also consider this effect by also excluding the 
p(1450) from the analysis. The p(1700) is a little too high for our region of interest. 

An understanding of the 7 = 0, 0++ channel has been elusive, despite much 
work. We have already included the low mass scalar (7(550) in order to understand 
TTTT scattering at low energy. In the 1 GeV region we expect the narrow resonance 
/o(980) to play a key role. The /o(980) has a 22% branching ratio into KK even 
though its central mass is below the KK threshold. The PDG also lists the scalar 
/o(1300) which has about a 93% branching ratio into tttt and a 7% branching ratio 
into KK. We shall use the /o(1300) here. It is hard to understand why, if the /o(980) 
is the ss member of a conventional 0+"'" nonet, it is lighter than the /o(1300). Most 
likely, the /o(980) is an exotic or a KK molecule [31]. If that is the case, its coupling 
to two pions ought to be suppressed in the 1/Nc picture. 

Now we will give, in turn, the tttt scattering amplitudes due to the exchange of 
the /o(1300), the /2(1270) and the p(1450). 

3.2.1 The Tensor /2(1270). 

In the first chapter we already observed that we represent the 3x3 matrix of tensor 
fields by a symmetric traceless matrix T^i,. T^^, transforms covariantly under a generic 
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chiral transformation (Eq. (1.26)). The 2-pion coupling can be deduced from the 
following chirally preserving term obtained from Eq. (1.27) 

-^{f2),Ad'^-d''n] . (3.6) 

In this case we note that the chiral invariant interaction is just the same as the minimal 
one we would have written down without using chiral symmetry. The partial width 
is then 

r(/.(1270) ^ .X) = , (3.7) 

where is the pion momentum in the /2 rest frame. This leads to I72I = 13.1 GeV~^. 

To calculate the /2 exchange diagram we need the spin 2 propagator [44] (see also 
Eq. (A.16)). 

?: n 11 

(3.8) 



g2 - ml 



2 (^Wl^M2t'2 ~l~ ^/ilf2^M2i^l) 2^/Ul/i2^!^lt'2 



where 



V = + ^ • (3-9) 

A straightforward computation then yields the /2 contribution to the tttt scattering 
amplitude: 

. / . N 7i f 16 4 10 2 1 2 I/.2 , 2N 



2 mts^ s*^ 



3 m^f Gmi 

J2 J2 J2 , 



(3.10) 



We notice that in the previous expression the behavior at high energy is dominated by 
the power. This is due to the fact that the present tensor is a massive particle not 
protected by any symmetry. This indicates that at higher energy we need to include a 
larger number of resonant states in order to bring the amplitudes inside the unitarity 
bounds. From this point of view our model is very close to the string mode [11], where 
an infinite number of resonances are identified with the string vibrational modes. The 
singularity in the propagator will be regulated as prescribed in the first chapter. 
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3.2.2 The vector meson yo(1450) and the /o(1300) scalar. 

The contribution of the s— wave radially excited p(1450) meson to the scattering 
amplitude is 



,2 



Ap>{s,t,u)- 



t{u — s) u{t — s 



+ 



m?pi — t m^, — u 



where Qp't^t^ is related to the p(1450) — > tttt partial width by 



(3.11) 



r(p(1450) - tttt) = 1^ . (3.12) 

With a branching ratio of 100% into two pions we get Igp'Tnrl — 7.9. In order to get 
the invariant amplitude in Eq. (3.11) there is no need to include the p(1450) as a 
massive chiral gauge field. [24] . 

By using Eq. (1.25)) we can easily deduce the chiral coupling to two pions 

+ ^f,d,7r-d^7r . (3.13) 



The partial decay width is 



r(/„(1300) ^ = 1 - 5^ X (M^ - 2m;)^ . (3.14) 

Since these resonance parameters are not well defined, for definiteness we will as- 
sume the PDG central values [37], i.e. rtot(/o(1300)) = 0.275 GeV and M/„ = 1.3 
GeV. Hence we deduce |7o| — 2.88 GeV~^. The unregularized invariant scattering 
amplitude is 

^/o(«,^,«) = |^^^^- (3.15) 
We will regularize this resonance propagator as for the /o(980) 



Ami 



3.2.3 /2(1270) + /o(1300) + p(1450) 

Now we are in a position to appraise the contribution to Rq of the next group of reso- 
nances. In order to better understand the local cancellation mechanism in this energy 
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range we will not consider the background effect at the moment. The contributions of 
each resonance are shown in Fig. 3.2 Note that the /o(1300) piece is not the largest, 

1 .5 

1 .0 

0.5 
°q;° 0.0 
-0.5 



-1 .0 
-1 .5 

0.3 0.8 1.3 1.8 

x/s GeV 

Figure 3.2: Contributions for /?[]. Solid line: f2(t + u). Dashed line: /2(s). 
Dotted hue: /o(1300). Dot-dashed line: p(1450) 

as one might at first expect. That honor goes to the /2 contribution which is shown 
divided into the s-channel pole piece and the {t + u) pole piece. We observe that the 
s-channel pole piece, associated with the /2, vanishes at = Mf^. This happens 
because the numerator of the propagator in (A. 16) is precisely a spin-2 projection 
operator at that point. The p(1450) contribution is solely due to the t and li— channel 
poles. It tends to cancel the t and channel pole contributions of the /2(1270) but 
does not quite succeed. The t and -u— channel pole contributions of the /o(1300) turn 
out to be negligible. Notice the difference in characteristic shapes between the s and 
{t + u) exchange curves. Fig. 3.3 shows the sum of all these individual contributions. 
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There does seem to be cancellation. At the high end, Rq starts to run negative well 
past the unitarity bound around 1.5 GeV. But it is reasonable to expect resonances 
in the 1.5 — 2.0 GeV region to modify this. The maximum positive value of Rq is 
about 1 at ^/s = 1.2 GeV. This would be a problem if we did not have the low energy 
contributions {current algebra +p+a). The background contribution is providentially 
negative (Fig. 2.8), showing once again the need for the a. Let us now check that 
1 .5 

1 .0 

0.5 

0.0 
-0.5 
-1 .0 
-1 .5 

0.3 0.8 1.3 1.8 

Vs GeV 

Figure 3.3: Sum of the contributions in Fig. 3.2. 

the next group does not essentially modify the low energy results up to 1.2 GeV. The 
somewhat positive net contribution of these resonances to Rq is compensated for by 
readjustment of the parameters describing the low lying sigma. 

It may be interesting to include the effect of the background phase for the /o(1300) 
as we have just seen that it was very important for a proper understanding of the 
/o(980). To test this possibility we reversed the sign of the /o(1300) contribution and 




46 



CHAPTER 3. EXPLORING THE 1 GEV REGION 



show the result as the sohd curve in Fig. 3.4. This sign reversal is reasonable since 
our model suggests a background phase of about 270° in the vicinity of the /o(1300). 
It can be seen that there is now a significantly greater cancellation of the next group 
particles among themselves up to about 1.2 GeV. The resulting total fits are shown 

1 .5 
1 .0 
0.5 



0.0 
-0.5 

0.80 1.00 1.20 1.40 1.60 

\|T (GeV) 

Figure 3.4: Contribution from the next group of resonances; the solid line 
is obtained with the reverse sign of the /o(1300) piece; the dashed line is 
as in Fig. 3.3. 

in Fig. 3.5 for assumed /o(980) — > tttt branching ratios of both 100% and 78.1% and 
the parameters associated with the fits are shown in Table 3.2. It is clear that the 
fitted parameters and results up to about 1.2 GeV are very similar to the cases when 
the next group was absent. Above this region there is now, however, a positive bump 
in Rq at around 1.3 GeV. This could be pushed further up by choosing a higher 
mass (within the allowable experimental range) for the /o(1300). Resonances in the 
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0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Ys (GeV) Ys" (GeV) 

(a) (b) 

Figure 3.5: Prediction for R^ with the next group of resonances, (a) 

assumes (column 5 in Table 3.2) (5i?(/o(980) ^ 27r) = 100%) while (b) 

assumes (column 6) (Bi?(/o(980) ^ 2%) = 78.1%). 



1500 MeV region, which have not been taken into account here, would presumably 
also have an important effect in the region above 1.2 GeV. Clearly there is not much 
sense, at the present stage, in trying to produce a fit above 1.2 GeV. In the last 
column of Table 3.2 we have neglected the effect of the p(1450). The resulting fit is 
shown in the last column of Table 3.2 and it is seen to leave the other parameters 
essentially unchanged. 

It thus seems that the results are consistent with the hypothesis of local cancel- 
lation, wherein the physics up to a certain energy E is described by including only 
those resonances up to slightly more than E and it is furthermore hypothesized that 
the individual particles cancel in such a way that unitarity is maintained. 
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3.3 Inelastic effects 



Up to now we have completely neglected the effects of coupled inelastic channels. Of 
course the An channel opens at 540 MeV, the Gn channel opens at 810 MeV and, 
probably most significantly, the KK channel opens at 990 MeV. We have seen that 
a nice undestanding of the tttt elastic channel up to about 1.2 GeV can be gotten 
with complete disregard for inelastic effects. Nevertheless it is interesting to see how 
our results would change if experimental data on the elasticity parameter rj^ are 
folded into the analysis. Figure 3.6 illustrates the results for r^o(*) obtained from an 

1 .2 n 1 1 1 1 1 1 1 1 1 1 ^ 



Tlo 



1 .0 



0.6 



0.4 



0.2 



0.0 



0.99 1.09 1.19 1.29 1.39 1.49 1.59 

\T (GeV) 



Figure 3.6: An experimental determination of t^q = wl — 4|T{2oP [45] 



experimental analysis [45] of tttt KK scattering. For simplicity, we approximated 
the data by a constant value rj^ = 0.8 above the KK threshold. Figure 3.7(a) shows 
the effect of this choice on Rq{s) computed without the inclusion of the next group 
of resonances, while Fig. 3.7(b) shows the effect when the next group is included. 
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0.2 0.4 0,6 0.8 1.0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 



Ys (GeV) Ys 
(a) (b) 

Figure 3.7: Predictions with phenomenological treatment of inelasticity 

{tJq = 0.8) above KK threshold, (a): without next group, (b): with neod 

group. 

Comparing with Fig. 3.1(b) and 3.5(b), we see that setting tjq — 0.8 has not made 

any substantial change. The parameters of the fit are shown in Table 3.2 as are the 
parameters for an alternative fit with rj^ = 0.6. The latter choice leads to a worse fit 
for 

We conclude that inelastic effects are not very important for understanding the 
main features of tttt scattering up to about 1.2 GeV. However, we will discuss the 
calculation of r)Q{s) from our model in the tttt — > KK paragraph. 

3.4 Phase Shifts 

Strictly speaking the Chiral Resonance Model only entitles us to compare the real 
part of the predicted amplitude with the real part of the amplitude deduced from 
experiment. Since the predicted Ro{s) up to 1.2 GeV satisfies the unitarity bound 
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(within the fitting error) we can calculate the imaginary part /o(s), and hence the 
phase shift Sq{s) on the assumption that full unitarity holds. This is implemented 
by substituting Rq{s) into Eq. (2.1) and resolving the discrete sign ambiguities by 
demanding that Sq{s) be continuous and monotonically increasing (to agree with 
experiment). It is also necessary to know f]Q{s) for this purpose; wc will be content 
with the approximations above which seem sufficient for understanding the main 
features of tttt— scattering up to 1.2 GeV. 

In this procedure there is a practical subtlety already discussed at the end of 
section IV of Ref. [8]. In order for S^{s) to increase monotonically it is necessary 
that the sign in front of the square root in Eq. (2.1) change. This can lead to a 
discontinuity unless 2|i?Q(s)| precisely reaches tIq{s). However the phase shift is rather 
sensitive to small deviations from this exact matching. Since the fitting procedure 
does not enforce that |-Ro(*)l precisely to 77o(s)/2 0.5, this results in some small 
discontinuities. (These could be avoided by trying to fit the phase shift directly.) 

Figure 3.8 shows the phase shift 5q{s) estimated in this manner for parameters 
in the first column of Table 3.2. As expected, the agreement is reasonable. A very 
similar estimate is obtained when (column 3 of Table 3.2) ?7q is taken to be 0.8 
while considering the tttt branching ratio of /o(980) to be its experimental value of 
78.1%. It appears that these two parameter changes are compensating for one other 
so that one may again conclude that the turning on of the KK channel really does 
not have a major effect. When the next group of resonances is included (column 7 
of Table 3.2) the estimated Sq{s) is very similar up to about 1.2 GeV. Beyond this 
point it is actually somewhat worse, as we would expect by comparing Fig. 3.7(b) 
with Fig. 3.7(a). 
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fV (GeV) 

Figure 3.8: Estimated phase shift using the predicted real part and uni- 
tarity relation. 
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3.5 The inelastic channel nn KK 

We have seen that tttt — > tttt scattering can be understood up to about 1.2 GeV 
without including this inelastic channel. In particular, a phenomenological description 
of the inelasticity did not change the overall picture. However we would like to begin 
to explore the predictions of the present model for this channel also. The whole 
coupled channel problem is a very complicated one so we will be satisfied here to 
check that the procedure followed for the tttt elastic channel can lead to an inelastic 
amplitude which also satisfies the unitarity bounds. Specifically we will confine our 
attention to the real part of the I — J — tttt ^ KK amplitude, i??2 o defined in 



In exact analogy to the tttt — > tttt case we first consider the contribution of the 
contact plus the K*{S92) plus the (7(550) terms. It is necessary to know the couphng 
strength of the a to KK, defined by the effective Lagrangian piece 



If the a is ideally mixed and there is no OZI rule-violating piece we would have 
'y^K'K ~ 7o as defined in Eq. (1.25). For definiteness, we shall adopt this standard 
mixing assumption. The appropriate amplitudes are listed in Appendix A. Figure 3.9 
shows the plots of i??2;o current algebra part alone, the current algebra plus 

K* and the current algebra plus K* plus a parts. Notice that unitarity requires 



The current algebra result already clearly violates this bound at 1.05 GeV . As before, 
this is improved by the K* vector-meson exchange contribution and further improved 
by the very important tail of the a contribution. The sum of all three shows a structure 



Eq. (A.ll). 



1(tKK 

2 



ad^Kd^K . 



(3.16) 




(3.17) 
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NTs' (GeV) 

Figure 3.9: Contributions to tttt KK (i??2;o)- The solid line shows the 
current algebra result, the dashed line represents the inclusion of K*{S92), 
the dotted hne includes the (t(550) too. 
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similar to the corresponding Fig. 2.8 in the tttt — > tttt case. The unitarity bound is 
not violated until about 1.55 GeV. 

Next, let us consider the contribution of the /o(980) which, since the resonance 
straddles the threshold, is expected to be important. We need to know the effective 
coupling constant of the /o to tttt and to KK. As we saw in Eq. (3.3), and the 
subsequent discussion, the effective tttt coupling should be taken as ^f^^^e^^. Ex- 
perimentally, only the branching ratios for /o(980) tttt and /o(980) — > KK are 
accurately known. We will adopt for definiteness the value of 7/o7r7r corresponding to 
the fit in the third column of Table 3.2 (rtoi(/o(980)) = 76 MeV). It is more difficult 
to estimate the /o(980) —>■ KK effective coupling constant since the central value of 
the resonance may actually lie below the threshold. By taking account ^ of the finite 
width of the /o(980) we get the rough estimate {'yf^xKl — 10 GeV~^ 4|7/(,7r7r| for 
the choice in the third column, M/q(98o) = 990 MeV. Of course, this estimate is very 
sensitive to the exact value used for Mjq(98o). It seems reasonable to take 7/„a'1? 
purely real. The results of including the /o(980) contribution, for both sign choices 
of ^f^KK-i shown in Fig. 3.10. The unitarity bounds are satisfied for the positive 
sign of 'Jf^^K-K but slightly violated for the negative sign choice. 

Finally, let us consider the contributions to tttt KK from the members of the 

^With rtot(/o(980)) = 76 MeV we would have r(/o(980) ^ KK)) = 16.6 MeV. Then 7^^^;^ is 
estimated from the formula: 



where A{f(){M) KK) is the reduced amplitude for an /o of mass M to decay to KK, $(M) is 
the phase space factor and p{M) is the weighting function given by 



16.6 MeV = 





Here, Mq is the central mass value of the /o(980). 
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Figure 3.10: Effect of /o(980) on tttt KK. The solid curve corresponds 
to a negative Jf^xK the dashed one to a positive sign. 
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multiplets containing the next group of particles. There will be a crossed channel 
contribution from the strange excited vector meson ir*(1410). However it will be 
very small since ir*(1410) predominately couples to K*7r and has only a 7% branching 
ratio to Kn. In addition there will be a crossed channel scalar _K'q(1430) diagram as 
well as a direct channel scalar /o(1300) diagram contributing to nn KK. The 
/o(1300) piece is small because /o(1300) has a very small branching ratio to KK. 
Furthermore the irQ(1430) piece turns out also to be small; we have seen that the 
crossed channel scalar gave a negligible contribution to tttt — > tttt. The dominant 
next group diagrams involve the tensor mesons. Near threshold, the crossed channel 
iir2(1430) diagram is the essential one since the direct channel /2(1270) contribution 
for the J = partial wave is suppressed by a spin-2 projection operator. Above 
1270 MeV the /2(1270) contribution becomes increasingly important although it has 
the opposite sign to the crossed channel tensor piece. Figure 3.11 shows the net 
prediction for -R?2 obtained with the inclusion of the main next group contributions 
from the ^"2(1430) and /2(1270). Both assumed signs for ^f^xK shown and 
other parameters correspond to column 3 of Table 3.2. Clearly there is an appreciable 
effect. Figure 3.12 shows the magnitude of |-R?2;ol together with one experimental 
determination [45] of |T{'2.oI = \/(-Ri2;o)^ + (-^i2;o)^- The positive sign of 'Jf^^xK 
favored but, considering the uncertainty in |7yjj;^xl among other things, we shall not 
insist on this. It seems to us that the main conclusion is that the unitarity bound 
can be satisfied in the energy range of interest. In this analysis we have shown that 
the local cancellation principle is satisfied and that the a plays an important role in 
the TTTT — > KK channel. 
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Figure 3.11: Effects on tttt KK due to tfie next group of resonances for 
the two different sign clioices in Fig. 3.10. 
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3.6. THE LONELY a 

3.6 The lonely a 
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In Reference [33], Tornqvist and Roos presented a model of tttt scattering which 
supports the existence of the old a meson at a pole position, Sq^^ = 0.470—^0.250 GeV. 
While this model is constructed to satisfy unitarity, it does not explicitly take crossing 
symmetry into account. In particular, one may question [34] the validity of neglecting 
the crossed-channcl p meson exchange contributions, which arc generally considered to 
be important. It is actually very complicated, as noted by the authors themselves, to 
examine this question in their model. We can investigate this issue in the framework 
of the Chiral Resonance Model [8-10, 55] We find that the consistent neglect of the 
p exchange docs not destroy the existence of the a meson but merely modifies its 
parameters so that they get close to the results of Ref . [33] . 

In the previous chapters [10] a best fit to the real part of the I — J — partial 
amphtude was found for a mass — 559 MeV, a width G' — 370 MeV and 

1 /2 

G/G' = 0.29 (pole position Sq = 0.585 — i0.176 GeV). It is an easy matter to 
neglect the p meson contributions (including the associated contact term needed for 
chiral symmetry) and make a new fit. The resulting Rq in comparison with the 
experimental data is shown in Fig. 3.13 and is about as good as the previous fit 
including the p meson. (Of course, the p meson is definitely present in nature.) The 
new fitted parameters are the mass M„ = 378 MeV, the width G' = 836 MeV and 

1 /2 

G/G' = 0.08 (pole position Sq = 0.495 — i0.319 GeV). The new mass and width are 
close to the values found in Ref. [33] . We therefore would expect that including the 
p exchange in their framework would raise their mass by roughly 100 MeV and lower 
their width prediction. This behavior can be easily understood in a qualitative sense, 
since the addition of the p raises the energy at which the unitarity bound is violated 
(see Fig. 2.6). Of course, in the previous chapters, the question of whether the a and 
/o(980) are gg, q^q^ states or some superposition is not directly addressed. 
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3.7. CONCLUSIONS FOR THE CHIRAL RESONANCE MODEL 61 

3.7 Conclusions for the Chiral Resonance Model 



We have used the Chiral Resonance scheme to obtain a simple approximate analytic 
form for the real part of the tttt— scattering amplitude in the energy range from 
threshold to about 1.2 GeV. It satisfies both crossing symmetry and (more non- 
trivially) unitarity in this range. Inspired by the leading 1/Nc approximation, we 
have written the amplitude as the sum of a contact term and poles. Of course the 
leading 1/Nc amplitude can not be directly compared with experiment since it is 
purely real (away from the direct channel poles) and diverges at the pole positions. 
Furthermore, an infinite number of poles and higher derivative interactions are in 
principle needed. To overcome these problems we have employed the following Chiral 
Resonance Model procedure: 

a. We specialized to predicting the real part of the amplitude. 

b. We postulated that including only resonances from threshold to slightly more 
than the maximum energy of interest is sufficient. We have seen that this local 
cancellation appears stable under the addition of resonances in the 1300 MeV 
range. Beyond this range we would expect still higher resonances to add in such 
a way as to enforce unitarity at still higher energies. 

c. In the effective interaction Lagrangian we included only terms with the minimal 
number of derivatives consistent with the assumed chiral symmetry. 

d. The most subtle aspect concerns the method for regularizing the divergences at 
the direct channel resonance poles. In the simplest case of a single resonance 
dominating a particular channel (e.g. the p meson) it is sufficient to add the 
standard width term to the denominator (e.g. the real part of Eq. (1.10)). 
For an extremely broad resonance (hke a needed low-energy scalar isosinglet) 
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the concept of width is not so clear and we employed the slight modification 
of the Breit-Wigner amplitude given in Eq. (1.11). Finally, for a relatively 
narrow resonance in the presence of a non-negligible background, we employed 
the regularization given in Eq. (1.12) which includes the background phase. 
Sclf-consistcncy is assured by requiring that the background phase should be 
predicted by the model itself. 

All the regularizations introduced above are formally of higher than leading order 
in the 1/Nc expansion (i.e. of order and higher) and correspond physically to 

pion rescattering effects. This rescattering is schematically represented in Fig. 3.14. 
In the case of non-negligible background phase, there is an interesting difference 

CK 

Figure 3.14: Effective resonance— 27r couphng due to the pions' rescatter- 
ing effects. The latter has been shown schematically in a generic pertu- 
bative scheme. For simplicity we only considered the rescattering due to 
a four pion contact term (•). 

from the usual tree-level treatment of pole diagrams. The effective squared couphng 
constant, gj^^^ of such a resonance to two pions, is then not necessarily real and 
positive. Since this regularization is interpreted as a rescattering effect it does not 
mean that ghost fields are present in the theory. This formulation maintains crossing 
symmetry which is typically lost when a unitarization method is employed. 
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In this analysis, the most non-trivial point is the satisfaction of the unitarity 
bound for the predicted real part of the partial wave amlitude, 

\Ri\<^, (3.18) 

where -qf < 1 is the elasticity parameter. The well-known difficulty concerns Rq. If 
ry/(s) is known or calculated, the imaginary part //(s) can be obtained, up to discrete 
ambiguities, by Eq. (2.1). 

The picture of tttt scattering in the threshold to shghtly more than 1 GeV range 
which emerges from this model has four parts. Very near threshold the current al- 
gebra contact term approximates -Ro(s) very well. The imaginary part Iq{s), which 
is formally of order 1/N^ can be obtained from unitarity directly using Eq. (2.1) or, 
equivalently, by chiral perturbation theory. At somewhat higher energies the most 
prominent feature is the p meson pole in the I — J — 1 channel. The crossed chan- 
nel p exchange is also extremely important in taming the elastic unitarity violation 
associated with the current algebra contact term (Fig. 2.6). Even with the p present. 
Fig. 2.6 shows that unitarity is still violated, though much less drastically. This 
problem is overcome by introducing a low mass 550 MeV, extremely broad sigma 
meson. It also has another desirable feature: Rq{s) is boosted (see Fig. 2.9) closer 
to experiment in the 400 — 500 MeV range. The three parameters characterizing 
this particle are essentially the only unknowns in the model and were determined by 
making a best fit. In the 1 GeV region it seems clear that the /o(980) resonance, 
interacting with the predicted background in the Rams auer- Towns end eSect manner, 
dominates the structure of the / = J = phase shift. The inelasticity associated 
with the opening of the KK threshold has a relatively small effect. However we also 
presented a preliminary calculation which shows that the present approach satisfies 
the unitarity bounds in the inelastic tttt — > KK channel. 
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Other recent works [32,42,49,51,56,57] which approach the problem in different 
ways, also contain a low mass broad sigma. The question of whether the lighter scalar 
mesons are of qq type or meson-meson type has also been discussed [32,42,56]. In 
our model it is difficult to decide this issue. Of course, it is not a clean question from 
a field-theoretic standpoint. This question is important for understanding whether 
the contributions of such resonances are formally leading in the 1/Nc expansion. We 
are postponing the answer as well as the answer to how to derive the rescattering 
effects that were used to regularize the amplitude near the direct channel poles as 
higher order in 1/Nc corrections. Presumably, the rescattering effects could some 
day be calculated as loop corrections with a (very complicated) effective Wilsonian 
action. This would be a generalization of the chiral perturbation scheme of pions. 
Another aspect of the 1/Nc picture concerns the infinite number of resonances which 
are expected to contribute already at leading order. One may hope that the idea of 
local cancellation will help in the development of a simple picture at high energies 
which might get patched together with the present one. Is the simple high energy 
theory a kind of string model ? 

In this first part of the thesis we have shown [8-10,55] that it is possible to build 
a reasonable model for light meson interactions, which wc have called the Chiral Res- 
onance Model (ChRM). Using the previous scheme we have demonstrated that it is 
possible to understand tttt— scattering up to the 1 GeV region by shoehorning together 
poles and contact term contributions and employing a suitable regularization proce- 
dure. It seems likely that any crossing symmetric approximation will have a similar 
form. We will regard the ChRM as a leading order l/N^ mean field approximation 
for the Quantum Chromodynamics. 



Part II 
Heavy Systems 
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Chapter 4 



Heavy Baryons in the Bound State 
Approach 

4.1 Brief Introduction to the Heavy Physics 

In the second part of this thesis we will investigate the heavy baryon spectra. A 
heavy hadron schematically consists of a heavy quark Q of spin -Sheavy = 1/2 and a 
hght cloud. The latter describes the hght degrees of freedom with total light spin jV 
There also exist heavy hadrons which contain a higher number of heavy quarks, but 
we will not consider them here. The total hadron spin is obtained by adding together 
the heavy quark spin and the spin of the light cloud 

J ^ Shea.vy+jl ■ (4-1) 

A heavy baryon corresponds to an integer value of j;, while a heavy meson to half-odd 
integer values of j^. 

In the heavy quark mass limit (mg/A ^1) the heavy spin commutes with the 
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QCD hamiltonian Hq 

[Ho, 5heavy] = . (4.2) 

Hence the heavy spin is now a good quantum number leading to the heavy spin 
symmetry. Since in the heavy hmit the hamiltonian cannot depend on the mass of 
the heavy quark it is not possible to distinguish among heavy hadrons made with 
different heavy quarks. We will indicate this symmetry as heavy flavor symmetry. 
Heavy spin symmetry predicts that for given ji ^ and fixed parity, we have a 
degenerate doublet of total spin 

J = Ji^l- (4-3) 

while for ji = we have a heavy baryon with total spin J = 1/2 which we can 
identify as Aq, where Q indicates the heavy quark contained in the heavy hadron. 
Let us explicitly demonstrate that the splitting between the — 1" and — 0~ 
states of a Qq meson must vanish in the limit of infinite heavy quark mass. Since the 
action of "S^eavy oil ^ 0~ state produces a 1~ state, i.e. \Mi- >= 25'heavylM)- > we 
then have 

Ho\Mi- >= mi- 1 Ml- >= 2Sl^^^yHo\Mo- >= mo- |Mi- > , (4.4) 

implying that mi- — mo- — > as mg — > oo. Experimentally [37] we have 

^^^^1^ ^ 8% , ^ 1% , (4.5) 

mo ruB 

which indicates the goodness of the 1/Af expansion. We also notice that the heavy 
spin breaking is 0{1/M) as can be understood by looking at a one gluon exchange 
diagram. 
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4.2 Bound State Approach to the Heavy Baryon 
system 

There has been recent interest in studying heavy baryons (those with the valence 
quark structure qqQ) in the bound state picture [59, 60] together with heavy quark 
spin symmetry [61]. In this picture the heavy baryon is treated [62-67] as a heavy 
spin multiplet of mesons (Qq) bound in the background field of the nucleon (qqq), 
which in turn arises as a soliton configuration of light meson fields. 

A nice feature of this approach is that it permits, in principle, an exact expansion 
of the heavy baryon properties in simultaneous powers of 1/M and In the 

simplest treatments, the light part of the chiral Lagrangian is made from only pion 
fields. However it has been shown that the introduction of light vector mesons [64-66] 
substantially improves the accuracy of the model. This is also true for the soliton 
treatment of the nucleon itself [68-70]. Furthermore finite M corrections as well 
as finite Nc (nucleon recoil) corrections are also important. This has been recently 
demonstrated for the hypcrfinc splitting problem [71,90] and it will be explained in 
some detail in the next chapter. 

Since the bound state-soliton approach is somewhat involved it may be worthwhile 
to point out a couple of its advantages. In the first place, it is based on an effective 
chiral Lagrangian containing physical parameters which are in principle subject to 
direct experimental test. Secondly, the bound state approach models a characteristic 
feature of a confining theory. When the bound system is suitably "stretched" it does 
not separate into colored objects but into physical color singlet states. 
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4.2.1 Effective Lagrangian for the Heavy-Light system 

Here we will review the chiral Lagrangian for the low lying heavy mesons in the heavy 
limit. In the next chapter we will see how to include next to leading corrections in 
1/M. The model is based on a chiral Lagrangian with two parts, 

^ = -flight + -^heavy (4.6) 

The light part involves the chiral field U = = exp (2i0/F^), where (p is the 3x3 
matrix of standard pseudoscalars. Relevant vector and pseudovector combinations 
are Vfj, and defined in Eq. (1.17). In addition light vector mesons are included in a 
3x3 matrix field p^. The explicit form of flight is as in Ref. [90]. 

The heavy field H, which describes the heavy multiplet which contains the heavy 
pseudoscalar P and a heavy vector meson is 

H = 7o^S, (4.7) 

where in the leading order in 1/M the pseudoscalar fiuctuation field (P') and the 
vector meson fiuctuation (Q^) are connected with the heavy fields via 

P = e-'^^-^P' , = e-'^^-^g; . (4.8) 

is the super selected four velocity of the heavy meson. 
Under a heavy spin transformation S, H transforms as 

H^SH , (4.9) 

while under a chiral transformation H transforms as a matter field 



H HK^ , 



(4.10) 
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where K is defined in Eq. (1.15). The heavy spin and chiral preserving lagrangian 
takes the form [76] 



Aeavy/M = iV^Tv 



HD^'H 



dTr 



H Tr 

my 



Hj^j,F^''{p)H , (4.11) 



where Di^ = d^ + iagp^ - i(l - a)v^, and F^^{p) = - + ig [p^ , p^] . Further- 
more, mv is the hght vector meson mass while d ~ 0.53 and c ~ 1.6 are respectively 
the heavy meson-pion and magnetic type heavy meson-light vector meson coupling 
constants; a is a coupling constant whose value has not yet been firmly established. 



4.3 Mechanics of Baryon States 

Following the Callan-Klebanov idea [59], we first find the classical solution of the 
light meson action and then obtain the classical approximation to the wave function in 
which this baryon as soliton is bound to a heavy meson (yielding a heavy baryon). The 
hedgehog ansatz for the classical light baryon in the SU (2) case simply corresponds 
to 

r i 

^x) =exp -x-tF{\x\) , 
Pic = ^^|^| Qfca^fcG'(|a;|) , 

pl = a;,c = , (4.12) 

where p^c = (^^^c + ''""P^c) ^ is a coupling constant. The appropriate boundary 
conditions are 



F(0) = -71 , G(0) = 2 , cu'(O) = , 
F(oo) = G{oo) = uj{oo) = , 



(4.13) 
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which correspond to unit baryon number. The wave function for the heavy meson 
bound to the background Skyrmion field (4.12) is conveniently presented in the rest 
frame, = 0. In this frame 

\ 

(4.14) 



He 







with a, I, h representing respectively the isospin, light spin and heavy spin bivalent 
indices. The calculation simplifies if we deal with a radial wave function obtained 
after removing the factor x ■ t: 

ht, = ^{^-r)^i'<ii,,, (4.15) 



where «(|a?|) is a radial wave function, assumed to be very sharply peaked near |a?| = 
for large M (i.e. r^|M(r)p ~ S{r)). The heavy spinor Xh can be trivially factored out 

^/^di,h = ^diXh (4.16) 

in this expression as a manifestation of the heavy quark symmetry. We perform a 
partial wave analysis of the generalized "angular" wave function ip^f 

{g. 93-, r, ^) = E C:Ss-,,sYrUk, ks) . (4.17) 

Here Y^^ stands for the standard spherical harmonic representing orbital angular 
momentum r while C denotes the ordinary Clebsch-Gordan coefficients. ^d/(fc, fcs) 
represents a wave function in which the "light spin" and isospin (referring to the 
"light cloud" component of the heavy meson) are added vectorially to give 

1^ = /light + -Slight , (4.18) 

with eigenvalues — k{k + 1). The total hght "grand spin" 

g^r + K, (4.19) 
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is a good quantum number in the heavy hmit. It is also convenient to define the total 
grand spin operator 

G^g + 5heavy , (4-20) 

where Sheavy is the spin of the heavy quark within the heavy meson. Substituting 
the wave-function (4.15) into / d^xjChea.vy given in Eq. (4.11) yields the potential 
operator 

V = J dn ijj* {(T ■TAi + lA2}'ijj 

= J dOV'*{4Ai5iighf /light + 1A2}V' 



2Ai 



'k{k + 1) - ^ 



+ A2 , (4.21) 



where / dQ is the solid angle integration and Eq. (4.18) was used in the last step. In 
addition 

Ai = \dF'{0)-^G"{0), 

A2 = -^^(0)- (4-22) 

The A2 term is relatively small [65,66,90]. Both terms in Ai are positive with the 
second one (due to fight vectors) sfightly larger. There are just the two possibifities 
A; = and k = 1. It is seen that the A; = states, for any orbital angular momentum 
r, will be bound with binding energy 3Ai. The k = 1 states are unbound in this 
limit. The parity of the bound state wave function is 

Pb = i-iy , (4.23) 

which emerges as a product of (—1)*^ for in Eq. (4.17), —1 for the x ■ t factor in 
Eq. (4.15) and —1 due to the fact that the mesons (P, Q) bound to the soliton have 
negative parity. 
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4.4 Collective Quantization 

In the soliton approach, the particle states with definite rotational and fiavor quantum 
numbers appear after the so-called rotational collective modes are introduced and the 
theory is quantized. This is conveniently done by first finding the time independent 
parameters which leave the theory invariant, then those collective parameters are al- 
lowed to depend on time. The collective angle-type coordinate A{t) is introduced [77] 
as 

ax,t) = A{t)Ux)A\t), 
T-p{x,t) = A{t)T-p,{x)A-\t) , 

H{x,t) = A{t)H,{x) , (4.24) 

where a-nd are defined in Eq. (4.12) and H^ in Eqs. (4.14) and (4.15). For our 
purposes the important variable is the "angular-velocity" f2 defined by 

A^A^^-T-^l , (4.25) 

which measures the time dependence of the collective coordinates A(t). It should 
furthermore be mentioned that, due to the collective rotation, the vector meson field 
components which vanish classically (pg and uii) get induced: 

1 

Ui = — —(p{r)eijk^jfk and pg = --^ (Ci(r)lljfc + 6(0^ ' ^h) ■ (4.26) 

Substituting Eq. (4.24) and Eq. (4.26) into / d^x {jCugtt + >Cheavy) gives an additional 
contribution to the lagrangian of the general form 

= l^y'^' -X^-G, (4.27) 

in which a and x represent spatial integrals over the profiles in Eq. (4.12). The in- 
duced radial functions (fi{r), ^i(r) and ^2(^) are obtained from a variational approach 
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to [85]. For each bound state solution He, there will be a tower of states character- 
ized by a soliton angular momentum J^°' and the total isospin I satisfying I — J^°^ 
The soliton angular momentum is computed from this collective Lagrangian as 

= ^ , (4.28) 
while the total baryon angular momentum is the sum 

J = ^ + J-l + 5heavy ■ (4.29) 

The rotational collective Hamiltonian is obtained by performing the standard 
Legendre transform 



+ . (4.30) 



The moment of inertia is identified from the light soliton sector as = |(mA — 
uin) in terms of the nucleon and A masses. Equation (4.27) can be simplified by 
noting that the total angular momentum is given by J = J^°^ + G. Then we deduce 
the heavy baryon mass formula [59] 

Hcoii = ^(^A - ttin) [(1 - x)/(/ + 1) + XJ(J + 1) + • • •] , (4.31) 

where the eUipsis stands for the x(x — 1)G{G + 1) term which does not split the 
heavy baryon masses. In the heavy limit >Cheavy leads to x = 0. Thus we have the 
final results 

m(S^)-m(SQ) = 0, (4.32) 

2 

m{EQ)-m{AQ) = -{rriA - , (4.33) 

wherein the subscript Q denotes the heavy baryon which contains the heavy quark 
Q. It may be interesting to compare the experimental determination [37] for m(Ec) — 
m(Ac) ~ 170 MeV whith the theoretical prediction in Eq. (4.33) which provides 195 
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MeV. The result tti^Eq) — m(EQ) = is, of course, expected in the heavy quark spin 
symmetry hmit. However if we consider heavy spin violating term (i.e. x 7^ 0) we get 

^(^q) - ^(Sq) = X {rriA - tun) ■ (4.34) 

The next chapter will be devoted to the calculation of the hyperfine splitting param- 
eter X- 



Chapter 5 

Heavy Baryon Hyperfine Splitting 



5.1 Introduction 

We have already noticed that a compeUing feature of the heavy sohton approach 
is that it permits, in principle, an exact expansion of the heavy baryon properties 
in simultaneous powers of 1/M, l/N^. and, since it is based on a chiral Lagrangian, 
number of derivatives acting on the light components of the heavy system. In practice 
there are obstacles related to the large number of unknown parameters which must be 
introduced. Rather than treating the light soliton in a model with many derivatives 
of the light pseudoscalar fields it turns out to be much more efficient to use the 
light vector mesons. Based on a model [65] of the light vector interactions with 
the heavy multiplet, the leading order (in the 1/Nc and 1/M expansions) heavy 
baryon mass splittings have been discussed [66], obtaining satisfactory agreement 
with experiment. Actually the need for light vector mesons is not surprising since, in 
the soliton approach, they are necessary to explain, for example, the neutron-proton 
mass difference [68] and the nucleon axial singlet matrix element [69]. 

In the present chapter we focus our attention on the hyperfine splitting, which is 
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of subleading order both in 1/M and This is a more comphcated calculation 

and also involves using a cranking procedure [77] to obtain physical states which 
carry good spin and isospin quantum numbers. The first calculation of the heavy 
baryon hypcrfinc splitting in the perturbative bound state framework was carried out 
by Jenkins and Manohar [62] who got the formula 

where M* — M is the heavy vector-heavy pseudoscalar mass difference, d is the light 
pseudoscalar-heavy meson coupling constant and -F'(O) is the slope of the Skyrme 
"profile function" at the origin. This formula is obtained (see also section 5) by using 
the leading order in number of derivatives (zero) and leading order in 1/M heavy 
spin violation term. Therefore it is expected to provide the dominant contribution. 
Unfortunately, on evaluation, it is found to provide only a small portion of the ex- 
perimental S*-Ec masss difference. This naturally suggests the need for including 
additional higher order in derivative heavy spin violation terms. However, there are 
many possible terms with unknown coefficients so that the systematic perturbative 
approach is not very predictive. 

To overcome this problem we employ a relativistic Lagrangian model [65] which 
uses ordinary heavy pseudoscalar and vector fields rather than the heavy "fluctuation" 
field multiplet [61]. This model reduces to the heavy multiplet approach in leading 
order and does not contain any new parameters. We will show that [71] such a 
model (considered, for simplicity, to contain only light pseudoscalars; i.e., the light 
part is the original Skyrme model [80]) yields a "hidden" heavy spin violation which 
is not manifest from the form of the Lagrangian itself. This hidden part involves 
two derivatives and is actually more important numerically than the zero derivative 
"manifest" piece which leads to Eq (5.1). However this new result is still not sufficient 
to bring the predicted S*-Ec mass difference into agreement with experiment. The 
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prediction for this difference is actually correlated to those for Ec-Ac and A-A^, the 
A - nucleon mass difference by [59]: 

m (E*) - m (Ee) =m{A)-m{N)-^ [m (EJ - m (A)] . (5.2) 

This formula depends only on the collective quantization procedure being used rather 
than the detailed structure of the model. If m (Ec) — m (A) and m (A) — m (N) are 
taken to agree with experiment, Eq (5.2) predicts 41 MeV rather than the experimen- 
tal value of 66 MeV. This means that it is impossible to exactly predict, in models of 
the present type, all three mass differences which appear in Eq (5.2). The goodness 
of the overall fit must be judged by comparing all three quantities with experiment. 
Our focus, of course, is the left hand side of Eq (5.2) which is of order 1/M while the 
right hand side involves the difference of two order M° quantities. A similar calcula- 
tion in the model with only light pseudoscalars was carried out by Oh and Park [67] . 
However, they did not make a 1/M expansion in order to reveal the hidden violation 
terms. They also introduced a one-derivative "manifest" heavy spin violation term 
with a new relatively large unknown constant in order to improve the agreement with 
experiment. 

In the present chapter we show that it is not necessary to introduce any new vio- 
lation terms to agree with experiment if a chiral Lagrangian including light vectors is 
employed. Typical results are summarized, compared with experiment and compared 
with the Skyrme model for the light sector in Table 5.1. A much more detailed discus- 
sion is given later in the text. We notice from the last row, that the model with hght 
vectors gives a very satisfactory account of the E*-Ec hyperfine splitting in contrast 
to the model without hght vectors. There are also noticeable effects when the use 
of the heavy meson reduced mass is taken as a simple approximation for kinematical 
corrections. Similarly, the first four rows of Table 5.1 show that the other predictions 
of the model with light vectors agree well with experiment. 



80 



CHAPTER 5. HEAVY BARYON HYPERFINE SPLITTING 



mass difference 


expt. 


present model 


present model + CM 


Skyrme 


Ac-N 


1345 


1257 


1356 


1553 


Ab - Ac 


3356±50 


3164 


3285 


3215 


A:.-Ae 


308 


249 


342 


208 


Se-A, 


168 


172 


158 


185 




66 


42 


63 


16 



Table 5.1: Typical results for the present model (including light vectors) 
compared with model with hght pseudoscalars only ("Skyrme" column) 
and compared with experiment. No "manifest" heavy spin violation ef- 
fects other than M* ^ M have been included. The column "present 
model + CM" simply takes into account recoil corrections by replacing 
the heavy meson mass by the reduced mass. A^ denotes a negative parity, 
spin 1/2 state. The quantity a in Eqs (5.4) was taken to be zero. All 
masses in MeV. 
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5.1.1 Relativistic Lagrangian for the Heavy Mesons. 

For the sector of the model describing the hght pseudoscalar and vector mesons we 
adopt the chirally invariant Lagrangian discussed in detail in the literature [86,87,90]. 
We now present the relativistic Lagrangian, which describes the coupling between the 
light and heavy mesons [65] 

+2iMd [Pp^Q^ - Q^p>^P) + "^e"^^^ [QuaP^Qp + QpP^Q^a] (5.3) 
+ ^ \2Q,F'^^ ip) - -e^^^'' [DpPF,, (p) + Q^F^, (p) DpP] | . 

Here we have allowed the mass M of the heavy pseudoscalar meson P to differ from 
the mass M* of the heavy vector meson Q^. Note that the heavy meson fields are 
conventionally defined as row vectors in isospin space. The covariant derivative in- 
troduces the additional parameter a: 

D, (P, g J = {d^ + tagp, -t{l-a)v,) (P, g J . (5.4) 

The covariant field tensor of the heavy vector meson is then defined as 

Qm- = D^Q. - DJQ^ . (5.5) 

The coupling constants d, c and a, which appear in the Lagrangian (5.3), have still 
not been very accurately determined. In particular there is no direct experimental 
evidence for the value of a, which would be unity if a possible definition of light vector 
meson dominance for the electromagnetic form factors of the heavy mesons were to 
be adopted [88]. We will later adjust a to the spectrum of the heavy baryons. The 
other parameters in (5.3) will be taken to be: 



d = 0.53 , c= 1.60 ; 
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M = 1865MeV , M* = 2007 MeV , 



D — meson ; 



M = 5279MeV , M* = 5325 MeV , 



B — meson. 



(5.6) 



It should be stressed that the assumption of infinitely large masses for the heavy 
mesons has not been made in (5.3). However, a model Lagrangian which was only 
required to exhibit the Lorentz and chiral invariances would be more general than 
the relativistic Lagrangian (5.3). Actually the coefficients of the various Lorentz 
and chirally invariant pieces in the relativistic Lagrangian (5.3) have precisely been 
arranged to yield the spin-fiavor symmetric model (4.11) in the heavy quark limit [76]. 

5.2 An apparent puzzle 

In this paragraph we will resolve an apparent puzzle which arises when calculating 
the corrections to the hyperfine splitting using the relativistic lagrangian presented 
in Eq. (5.3). Here we will neglect the light vector contributions which is equivalent 
to set a = and c = in Eq. (5.3). Then the heavy Lagrangian in Eq. (5.3) becomes 



Here we have modified the coefficient of the fifth term in Eq. (5.3) to include a new 
source of heavy spin breaking. First let us consider the calculation of the hyperfine 
splitting in the heavy field approach. This, of course, arises at first sub-leading order 
in 1/M and violates the heavy spin symmetry. Thus we must add to Eq. (4.11) 
suitable heavy spin violating terms [62] : 



+D^PD"F M^PP - -Q^uQ^" + M*^Q^Q' 
+ 2iMd (Pp^Q^ - Q^p^^P) 




(5.7) 



M - M* 



8 



Tr Ha^.Ha'"' + 



{d - d') 
2 



TV Hpi'H^^-f^ +••• . 



(5.8) 
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The first term has no derivatives while the second term has one derivative. The 
hyperfine sphtting is related to a collective Lagrangian parameter (see section 5.3 for 
details) x with a proportionality factor of the A-N mass difference: 

m(E^) - m(EQ) = [m(A) - m{N)] x ■ (5.9) 

(At present only Eg is well established experimentally.) For Eq. (5.8) we have 

M*-M d-d' 

The first term was obtained in Ref. [62] while the second seems to be new. Notice 
that (M* — M) and {d — d') behave as 1/M. These quantities are the same as the 
ones appearing in the ordinary field Lagrangian (5.7). It would thus seem that CJ^eaivy 
in Eq. (5.8) neatly summarizes the heavy spin violation in Eq. (5.7). 

Now let us consider the calculation of x from Eq. (5.7) directly based on exact nu- 
merical solution of the associated coupled differential equations. We content ourselves 
with the graphical presentation of some results^ and relegate the details to Ref. [90]^. 
Figure 5.1 shows x plotted against M for three cases: i) M* = M, d' — d — 0.53, 
ii) M* - M c:^ (0.258GeV)VM (a fit to experiment), d' = d = 0.53, iii) M* = M, 
d' — d — (0.0991GeV)/M (an arbitrary choice which sets the coupling constant split- 
ting to be 10% at the D meson mass). We immediately notice that x does not vanish 
when there is no manifest heavy spin violation, i.e., M — M*, d — d'. In fact the 
dominant part of the contribution to x for realistic heavy meson masses is already 
present in this case. By subtracting out this piece we note that the signs of the con- 
tributions due to M* ^ M and d' ^ d agree with those predicted in Eq. (5.10). It is 
interesting to note that all three curves in Fig. 5.1 fall off as 1/M for M > lOGeV. 
^For the Skyrme model parameters we use the experimental value of F„ and esk = 6.0. This 

results in a profile with F'(0) = 1.20 GeV. 

^Similar calculations were done in Ref. [67] but they did not consider the M = M*, d = d! case. 
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Figure 5.1: X vs. M computed by numerical integration. Solid line M* — 
M, d' = d; dotted line M* ^ M, d' ^ d, dashed line M* ^M,d' ^ d. 
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But our main task is to understand the source of the puzzhng non-zero contribution 
in case i. It is clear that the ordinary field Lagrangian (5.7) must contain heavy 
spin violating pieces which are not manifest. We will now explore this in detail by 
rewriting Eq. (5.7) in terms of the "fiuctuation field" H and expanding it in powers 
of 1/M. 

5.2.1 Expansion of Lagrangian 

Since the effects of M 7^ M* and d ^ d! were taken into account in Eq. (5.10) it 
is sufficient to expand Eq. (5.7) with M* — M and d' — d. To describe the heavy 
particle moving with four- velocity V^, we introduce the factorization 

P = e-^^^-^P' , = e-'^^-^g^ . (5.11) 

P' is the pseudoscalar "fluctuation field". is not exactly the vector fluctuation 
field since V ■ Q is, not constrained to be zero. We therefore introduce the correct 
fiuctuation field Q'^ by 

Q^.-Q'i, + V,V-Q , (5.12) 

which shows that V ■ Q' = Substituting Eqs (5.11) and (5.12) into the Lagrangian 
(5.3) gives, in addition to the leading terms of order M, the presently interesting 
terms of order M°: 

C{P,Q) = {order M)- P' 0^7^+ Q'^D^Q^ -Q'^D.Di^Q" 

+ MV -QV-Q + iM {d^Q'^V -Q-V- QD^Q^^ 

+ 2iMd [p'V -pV -Q-V -QV ■ pF^ + • • • , (5.13) 

where the three dots stand for terms of order 1/M. In contrast to the massless fields 
P' and Q', V ■ Q is seen to have the large mass M. We thus integrate it out using 
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the equation of motion 



(5.14) 



Substituting Eq. (5.14) back into Eq. (5.13) gives 



>C(P,g) = {order M)-P'D^P^ + Q'^D^Q^'' -iQ'^F^''{v)Wu 

- 2d [P'V ■ pD^Q^ + D^Q'^'V ■ p'P) 

- ide"^'^^ {dM'^p.Q', - Q'^ppD^W.) 
-M^P'{V ■pfP^+--- , 



(5.15) 



where F^i,{v) — d^Vy — d^v^ — i[vi^,Vi,]. In order to extract the heavy spin violating 
pieces it is convenient to rewrite the order M° Lagrangian in terms of the heavy 
multiplet field H Eq. (4.7). After some algebraic calculation we find 



C{H) 



heavy ~t" 2 



c 

+ d 



i-Tr 
8 



iTr 
2 



HD'^H 
D^HY^5{V-p)H 



[Hn.lu] F^^{v)H 



- -Tr 

4 



7 • DHp^H^j^ 



- -Tr 

4 



-Tr 



7 • DH^^^^pi'H 



a^^D^H-f'^V-p-f^a'^'H + h.c. 



+ d' 



-Tr 



H {V ■ pf h] + -Tr [a^.Ha'"' {V ■ pf H 



+ 



(5.16) 



where £hcavy is given in Eq. (4.11) with c = a = 0. At this stage we see that 
Eq. (5.16) actually contains pieces which are not manifestly invariant under the heavy 
spin transformations H — > SH, H — > HSK These pieces involve two derivatives. 



5.3 Hyperfine splitting from the Hidden Terms 

We now sketch the computation of the portion of x in Eq. (5.9) which results from the 
"hidden" heavy spin violation in Eq. (5.3) that has been made explicit in Eq. (5.16). 
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For this purpose one needs to collectively quantize the Lagrangian as described in 
paragraph 4.4. To leading order in M, the "angular part" of the ground state wave 
function in Eq. (4.15) is [64,65] 

V'dl,^ = -^edi52h ■ (5.17) 



The specific value of the index h results from the choice — G — 1/2 where G is the 
"grand spin" . To next leading order in M the ground state wave function receives a 
heavy spin violating admixture of 



',(2) 1 



^dl,h — 



h2 . (5.18) 



Finally, the hyperfine splitting parameter x is recognized by expanding the collective 
Lagrangian [59], in powers of f2 and picking up the linear piece L^ow = {x/'^)^3 + ■ ■ ■■ 
Noting that the A-nucleon mass difference is given by the moment of inertia, which 
relates the angular velocity to the spin operator [77], this piece of the Lagrangian 
yields Eq. (5.9) after canonical quantization of the collective coordinates [59]. There 
are two types of contribution to x- The first type, from the heavy spin violating terms 
proportional to d in Eq. (5.16), corresponds to the evaluation of heavy spin violating 
operators in the ground state (5.17). The second type corresponds to the evaluation 
of heavy spin conserving operators in the ground state which includes an admixture 
of Eq. (5.18) due to the Tr [-ff,-f^HF^''{v)H\ term in Eq. (5.16). The net result for 
the "hidden" part of x is 

x^mi^-D- (5.19) 



4M V 2dJ 

This equation is expected to hold for large M. To this should be added the "manifest" 
part given in Eq. (5.10). 

It is important to compare Eq. (5.19) with the result for x obtained by the exact 
numerical solution for the model based on Eq. (5.7). This is gotten as an integral over 
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the properly normalized radial functions <l>(r), . . . , '^2{r) which appear in the P-wave 
solution of the bound state equation [66] : 



pe"* . (5.20) 



The spinor p labels the grand spin of the bound heavy meson. The choice = +1/2 
corresponds to p = (1,0)''^. The heavy limit bound state wave function in Eq. (5.17) 
corresponds to the special choice 

$(r) oc u{r) , *i(r) = -$(r) , *2(?^) = -2$(r) and *4(r) = . (5.21) 

The numerical solution to the bound state equations exactly exhibits these relations 
for M,M* ^ oo [66]. 

Equation (5.19) has an interesting d-dependence and vanishes a,t d — l/-\/2, which 
actually is not too far from the experimental value of this quantity. In Fig. 5.2 we 
compare the rf-dependence of the exact numerical calculation with the perturbative 
result of Eq. (5.19). It is seen that the large M perturbation approach works rea- 
sonably well and the gross structure of the hyperfine splitting is reproduced. For 
a detailed comparison of the two treatments it is important to note that for fixed 
M = M* the binding of the heavy meson increases with d. In particular this implies 
that the wave function is only reasonably localized for large enough d. As a strong 
localization is a basic feature of the perturbative approach it is easy to understand 
why this calculation does not yield the exact (numerical) result for small d. In fact, 
as d increases the agreement expectedly improves. However, upon further increase of 
d (at finite M,M*), the numerical solution to the bound state equations shows no- 
ticeable deviations from the heavy limit relations (5.21), which causes the moderate 
differences at larger d. 
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Figure 5.2: The d dependence of x for M = M* = 30 GeV and d ^ d! . 
Solid line is the exact numerical calculation. Dashed line is the large M 
perturbation formula given in Eq. (5.19). 
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We have solved [71] the apparent puzzle associated with the use of a model La- 
grangian containing ordinary fields for computing the hyperfine splitting parameter 
X by carefully expanding the Lagrangian in powers of 1/M. The key point was the 
need to preserve the constraint V ■ Q' — ior the heavy vector fiuctuation field. 

Of course, such a model Lagrangian (which has been used in many calculations) 
is not exactly QCD. Nevertheless it seems reasonable since it automatically has the 
correct relativistic kinematics and satisfies the heavy spin symmetry at leading order. 
We have seen (Eq. (5.16)) that at next order in 1/M, it predicts the coefficients of 
many terms which otherwise would be unspecified by heavy spin symmetry. 

It is amusing to note that these 1/M suppressed terms involve two derivatives 
and are actually more important for the computation of x than the zero derivative 
term in Eq. (5.10). This is readily understandable since the dynamical scale in this 
calculation is the binding energy, m{B) +m{N) — m{Ai,) ~ 620 MeV which is rather 
large for neglecting light vector mesons, higher derivatives etc. [See, for example, 
Ref. [10].] 

We are regarding the Lagrangian (5.3) with a = c = as an illustrative model 
rather than as a realistic one for comparison with experiment. As indicated earlier 
it seems necessary to include, in addition to finite M corrections, the effects of fight 
vector mesons as well as nucleon recoil. The discussion of x in this more complicated 
model and further details of the present calculation will be given in the next paragraph 
[90]. 
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5.4 Perturbative Approach and the Vector contri- 
bution. 

The perturbative approach can illuminate several aspects of the hyperfine splitting 
problem. This is due to the heavy quark symmetry which is naturally exploited by 
making an expansion in powers of 1/M using the heavy field formalism. Our starting 
Lagrangian (5.3) has been set up in such a way as to yield a heavy quark symmetric 
result as M ^ oo when M = M* is assumed, cf. Eq. (4.11). The perturbative 
1/M expansion is more general (presumably exact) but less predictive. Thus the 
1/M expansion provides a useful calibration in the large M limit. Since it deals with 
perturbation matrix elements it provides us with a convenient classification of the 
various sources of hyperfine splitting. The method is also advantageous in that it can 
be extended, without too much algebraic work, to different channels of interest. On 
the other hand, once the particular channels of interest are settled on, it is clearly 
more convenient to employ the exact numerical solution, which efficiently sums up a 
class of 1/M corrections. 

The leading order Lagrangian (4.11) can be supplemented by terms which man- 
ifestly break the heavy quark symmetry to leading order (M" with the present nor- 
malization) as follows: 




H 




+ — Lt,U^^^hf^^''{p)H 



Tr Hpi^HY' 



'75 



(5.22) 



Here the (M — M*) term measures the heavy spin violation due to the heavy pseu- 
doscalar - heavy vector mass difference. The {d — d!) term measures the heavy spin 
violation induced by choosing different coefficients for the fifth and sixth terms in 
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Eq. (5.3), while the (c — c') term corresponds to choosing different coefficients for 
the last and next-to last terms in Eq. (5.3). Finally the a term corresponds to the 
leading term obtained by using different values of a in Eqs (5.4) for P and Q. Note 
that (M - M*), {d- d'), (c - c') and a all behave as 1/M. 

In addition to the terms in Eq. (5.22), which manifestly break the heavy quark 
symmetry, there are, in fact, "hidden" violation terms contained in Eq. (5.3). The 
explicit expression for the hidden terms in the model without light vectors is given 
in Eq. (5.16) [71]. These were shown to exist (for the model without light vectors) 
in Ref. [71] and arise from performing a detailed 1/M expansion of the relativistic 
Lagrangian. In Reference [90] the numerical study has confirmed that this is also 
true when light vector mesons are included. In the last paragraph we have shown (c/. 
Fig. 5.2) [71]) that the dependence on d of the hyperfine splitting computed from these 
hidden terms using the perturbative approach generally matched the exact numerical 
calculation. Hence wc shall not explicitly isolate the extra hidden terms due to the 
addition of the light vectors but shall content ourselves with the numerical treatment 
given in Ref. [90]. 

Here we discuss the computation of x in some detail. We have already noticed 
that the dynamics of the model dictates that the bound-states occur for /c = 0, in 
which case ^d/(0, 0) — e^ij\f2. For reader's convenience we display again the relevant 
wave functions. The bound-state wave-function simply is 

V'd/,^(0,0,0,0)-^eda/.- (5.23) 

VOTT 

The k—\ unbound wave-function with no orbital excitation (r = 0) is 

V'<ii,^(l,^3,0,l) = -^Ui{\9z)Xh . (5.24) 
When violations of the heavy quark symmetry are included, g is no longer a good 
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quantum number. The grand spin, is a good quantum number 



G' — g + iSheavy ■ 

In the notation of Eqs (5.23) and (5.24) we have the grand spin eigenstates 

= g;, = 1/2) = ^ 

= = 1/2) = 



(5.25) 



Stt 

1 



(5.26) 
.(5.27) 



Note that in Eq. (5.26) the G'^ = +1/2 wave function is 52h since the index 2 corre- 
sponds to +1/2 for the anti-quark wave-function. The two states (5.26) and (5.27) 
differ with respect to their g and K labels. 

Now let us consider the potential for the bound-state wave-function in the presence 
of the first heavy quark symmetry violating term in Eq. (5.22). Substituting the G'- 
eigenstates -0^^^ and -0^^^ from Eqs (5.26) and (5.27) into Eq. (4.11) and the first term 
of Eq. (5.22) yields, after a spatial integration, the potential matrix in the t/j'^^^-t/j^'^^ 
space: 



V 



dF'{0) 



I 



3 
-1 



+ 



M-M* 



( ^\ 

V3 2 



(5.28) 



where F{r) is defined in Eq. (4.12) and F' — dF/dr. The first matrix shows that ^'^^^ 
is bound while ■0^^^ is unbound in the heavy spin limit. Since the second matrix gives 
mixing between ij)'^^'^ and if)^"^^ the latter must be included in the presence of effects 
which break the heavy quark symmetry. The diagonalized bound wave function is 
seen to be 

(5.29) 



8 rfF'(O) 

This is the proper wave-function to be "cranked" in order to generate the heavy spin 
violation. Using it in Eq. (4.24), which is then substituted into the a = limit of the 
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first term of Eq. (4.11), contributes a term in tlie collective Lagrangian 

1 M*-M , , 

— iZa where y = — ; — t-^ . (5.30) 

By using the Wigner-Eckart theorem we may express this for states of either Gg as the 
matrix element of the operator xO • G'. For convenience we have chosen to consider 
our wave-function as representing the conjugate particle in Eq. (4.15). Hence the 
matrix element of G' in this section differs by a minus sign from that of G defined in 
J — G + J^°\ with Jf°^ = (dL/dQi), which is the appropriate one when we form the 
total heavy baryon spin. Then the collective Lagrangian, Lcou may be written (see 
section 4.4) 

i^cou = l^^^^ -X^ G (5.31) 

which again leads to the Hamiltonian (4.31) and hence to the well known formula cf. 
Eq. (5.1) 

m(E^) - m(EQ) = [m(A) - m{N)] x ■ (5.32) 

The purpose in deriving this again was to explain the perturbative method and our 
notation. 

Next we shall give some new perturbative "manifest" contributions to x from 
Eq. (5.22). When all these terms are included the potential V in Eq. (5.28) is modified 
so that the properly diagonalized wave-function which replaces Eq. (5.29) becomes 

iP^^^ + eiP^^^ , (5.33) 

with 

(M -M*) + ^{d- d!) F'{0) + V^au;{0) + ^i'^'^ ^'f^ 

There are two types of contribution to x- The first type is analogous to Eq. (5.30) and 
arises when Eq. (5.33) is cranked and substituted into Eq. (4.11). The second type 
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is obtained by substituting the leading order wave function ip^^^ into the (c — c') and 
a terms in Eq. (5.22). The complete expression for x resulting from the "manifest" 
heavy spin violation is 



X 



a 
+ 3 



:C-C' 



8 -2^2^(^1(0) -6(0)) -4V^ ip"{0). (5.35) 

TTiy 



The quantities ^(0), ^(0) and (p"{0) are defined in Eq. (4.26). This formula may 
be useful for quickly estimating the effects of heavy spin violation in the coupling 
constants, which were not explicitly given in the previous discussion. Unfortunately 
there is no determination of the magnitude of these effects from the mesonic sector 
at present. In the previous paragraphs [71] the discussion of the "hidden" heavy 
contributions to x was given for the Lagrangian with only light pseudoscalars. 

The hyperfine splitting just discussed is for the ground state or P-wave heavy 
baryons. It is of some interest to briefly consider the negative parity heavy baryons 
with one unit of orbital excitation (S-wavc). In the heavy spin limit these bound 
states correspond to the r = 1 and A; = choice in Eq. (4.15): 

M^,93,l,0) = ^Yrxh. (5.36) 

The spin, Jught of the "light cloud" part of the heavy baryon is gotten by adding this 
g = I piece to the sohton spin J^°'. For the / = (which implies J^°' = 0) heavy 
baryons one finds Jught = 1 and the degenerate multiplet 

{A'q(1/2),A'q(3/2)} . (5.37) 

For the / = J^°' = I heavy baryons, J^ght can be either 0, I or 2 and we find the 
degenerate heavy multiplets 

S'q(1/2) , 
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{e^(1/2),E'q(3/2)} , 

{e'q(3/2),E'q(5/2)} . (5.38) 

In general, the situation is even more complicated and the subject will be fully in- 
vestigated in the next chapter. At present there are experimental candidates [37] for 
a negative parity spin 1/2 baryon at 2593.6 ± 1.0 MeV and a negative parity spin 
3/2 baryon A; at 2626.4 ± 0.9 MeV. 

Since experimental information is available, it is especially interesting to consider 
the sphtting between the two A'g states in Eq. (5.37). This sphtting stems from the 
violation of the heavy quark symmetry. For the Aq type states the total spin coincides 
with the grand spin G so that Eq. (5.37) may be alternatively considered a G = 1/2, 
G = 3/2 multiplet. Since the good quantum number is G, we may in general expect 
the hyperfine parameter x to depend on G. The collective Hamiltonian takes the 
form 

H^on = ^ . (5.39) 

On general grounds we see that for the case of the Aq's the collective Hamiltonian con- 
tribution to the hyperfine splitting will be suppressed. Setting J^°' = in Eq. (5.39) 
shows that the hyperfine splitting is of order (x^) or equivalently of order (1/M^). 
Unlike the ground state which involves only the G = 1/2 P-wave channel, there is 
another possibility for hyperfine splitting here. It is allowed for the G = 1/2 and 
G — 3/2 bound state energies to differ from each other. In the Lagrangian with only 
light pseudoscalars this does not happen and the Aq(1/2) — Aq(3/2) splitting is of 
order 1/M^. However when light vectors are added, there are "hidden" 1/M terms, 
which violate the heavy quark symmetry as e.g. 



i Tr 



<ja^,H^,F^\p)D''H + h.c. . (5.40) 



This term is likely to generate splitting for the multiplet (5.37) to order 1/M by 
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giving different binding energies to the G = 1/2 and G — 3/2 channels. It would be 
very interesting to investigate this in more detail. 

Finally, we add a remark concerning an amusing conceptual feature in the com- 
putation of hyperfine sphtting among the five S'^'s in Eq. (5.38). The total angular 
momentum of each state is given by 

G 

rsol 



J = J'°'+5+'Sheavy , (5.41) 

J light 

where wc arc now considering each operator to be acting on the wave-function rather 
than its complex conjugate. We have illustrated two different intermediate angular 
momenta which can alternatively be used to label the final state. If Jiight is used, we 
get the heavy-spin multiplets in Eq. (5.38). On the other hand, when the hyperfine 
splitting is turned on, the choice G is convenient, because it remains a good quantum 
number. According to the laws of quantum mechanics, we cannot simultaneously use 
both to specify the states, since the commutator 



•^Ught ' =4i Jiight -(SheavyXflr) (5.42) 



is generally non-vanishing. This means that we cannot uniquely trace the splitting 
of, say, the |Eq(1/2), Eq(3/2)| heavy multiplet in Eq (5.38), as hyperfine splitting 
interactions are turned on. Physically, this causes a mixing between the E'g's of the 
same spin. Rather, we must look at the whole pattern of the five masses. On the 
other hand, the problem simplifies for the computation of the ground state hyperfine 
splitting in Eq. (5.32). In that case the bound state wave function is characterized 
by flf = 0. Thus the commutator in Eq. (5.42) vanishes, and it is "trivially" possible 
to track the hyperfine splitting as a mass difference. 
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5.5 Numerical Results 



In this section we will briefly comment on the numerical results obtained in [90] for the 
masses of the heavy baryons within the relativistic Lagrangian model discussed above. 
The numerical procedure requires the solution of coupled inhomogeneous differential 
equations; the details are provied in Ref . [90] . In particular we will concentrate on the 
spin and isospin splitting in the realistic case of finite heavy meson masses (5.6). It 
should be noted that sizable quantum corrections occur for the classical soliton mass 
Mci [89]. It seems that these corrections are (approximately) equal for all baryons. 
Hence we will only consider mass differences between various baryons. In that case 
the absolute value of the classical mass M^i is redundant. The parameters in the light 
sector can completely be determined from properties of the corresponding mesons 
and by the |^ and |^ light baryons [90]. The corresponding moment of inertia is 
Q/2 — 5.00GeV~^. In the preceding paragraphs Ref. [71] we have shown (in the case 
without light vectors) that a major fraction of the P-wave hyperfine constant is due 
to terms in the relativistic Lagrangian (5.3) which do not manifestly break the heavy 
spin symmetry rather than to terms, which explicitly break this symmetry; as for 
example M ^ M*. The numerical results in Ref. [90] provide a quantitative extimate 
of the hidden contribution, when the light vectors are also included, by performing the 
calculation using identical masses from the charm sector i.e. M = M* = 1.865GeV 
and furthermore a — 0.3. All other parameters are as in Eq. (5.6). This results in 
Xp — 0.080 (x for the P-wave). Prom table 5.2 we recognize that this is about 80% 
of the value obtained using the physical masses M — 1.865GeV and M* — 2.007GeV. 
In the case of the S-wave the hidden piece is even more dominant. Por the symmetric 
choice of the mass parameters one finds xs = 0.175 (x for the S-waves) which is 
more than 90% of the value displayed in Table 5.2 for a — 0.3. It is also possible to 
show [90] that the light vector model predicts a substantially larger xp- 
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Since the contribution of the manifest (M* — M) breaking term is relatively small 
it is reasonable to expect that the others will be small too. 

Let us next discuss the spectrum of the baryons containing a single heavy quark. 
For this case we assume the reahstic masses as in Eq. (5.6). In Table 5.2 the numerical 
results for the lowest S- and P-wave bound states in the charm sector are displayed. 
As already noted in Ref. [65] the binding energy (a;) decreases with growing coupling 
constant a. This is the case for both the P- and S-wave channels. For M — > oo the 
heavy limit (see Eq. (4.22)) 

cu IdF'iO) + -^G"{0) + ^uj{0) = 3Ai - A2 (5.43) 
2 gniy v2 

will be attained. We note that the hyperfine parameters in these two channels behave 

oppositely as functions of a. Here we have chosen to measure the mass differences 

with respect to the lightest charmed baryon, Ac. Hence the mass differences with 

respect to Ec and E* directly reflect the a-dependence of hyperfine parameter xp 

while the corresponding dependence of the binding energy up can be extracted from 

the splitting relative to the nucleon. In addition the splitting with respect to the 

negative parity charmed baryons reflects the a-dependence of the S-wave channel 

binding energy cus- Finally the mass difference to A;, contains the energy eigenvalues 

and hyperfine parameters computed with the B and B* meson masses in Eq. (5.6). 

While the mass difference to the nucleon is improved with a positive value for a, 
the agreement for the mass differences between the heavy baryons slightly deteriorates 
when increasing this parameter. Nevertheless, fair agreement with the experimental 
data is achieved for quite a range of a. 

Table 5.2 also contains the model predictions when the background sohton is 
taken from the basic Skyrmc model [77, 80] which does not include the light vector 
mesons. Here we have adjusted the only free parameter (eskyrme — 4.25) to reproduce 
the A-nucleon mass difference. From the Ac-nucleon mass difference we observe 
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a 


-0.1 


0.0 


0.1 


0.2 


0.3 


Expt. 


Skyrme 


UJp 


564 


544 


522 


500 


478 




243 


Xp 


0.147 


0.140 


0.131 


0.123 


0.114 




0.053 


UJs 


316 


298 


281 


264 


247 




57 


Xs 


0.172 


0.181 


0.189 


0.197 


0.205 




0.346 




171 


172 


174 


175 


177 


168 


185 


K 


215 


214 


213 


212 


211 


233 


201 


K 


250 


249 


245 


242 


238 


308 


208 


K 


415 


413 


408 


402 


397 


? 


335 


K* 


468 


467 


464 


461 


458 


? 


437 


N 


-1237 


-1257 


-1278 


-1299 


-1321 


-1345 


-1553 




3160 


3164 


3167 


3170 


3173 


3356 ± 50 


3215 



Table 5.2: Parameters for heavy baryons and mass differences witfi respect 
to Ac. Primes indicate negative parity baryons, i.e. S-wave bound states. 
All energies are in MeV. 
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a 


-0.1 


0.0 


0.1 


0.2 


0.3 


UJp 


811 


786 


762 


737 


713 


Xp 


0.055 


0.053 


0.050 


0.048 


0.045 


UJs 


639 


617 


595 


573 


552 


Xs 


0.043 


0.046 


0.049 


0.052 


0.055 




189 


189 


190 


190 


191 




206 


205 


205 


205 


205 


K 


171 


168 


167 


164 


161 




363 


359 


358 


354 


351 




375 


373 


371 


369 


367 


N 


-4397 


-4422 


-4446 


-4471 


-4494 



Table 5.3: Parameters for heavy baryons and mass differences with respect 
to A;,. Primes indicate negative parity baryons, i.e. S-wave bound states. 
All energies arc in McV. The empirical value for the relative position of 
the nucleon is 4701 ± 50MeV [37]. 

that in comparison with the nucleon the masses of the heavy baryons are predicted 
about 200 MeV too large. This confirms the above statement that the spectra of 
both the light and the heavy baryons can only be reasonably reproduced when light 
vector mesons are included. This conclusion can already be drawn from the too small 
binding energies [66]. The hyperfine corrections make only minor changes in the 
Afo-Ac splitting. 

In Table 5.3 we display the analogous predictions for the bottom sector. According 
to the heavy spin symmetry the binding energies of the P- and S-wave channels 
approach each other. Hence the mass differences between the even and odd parity 
baryons containing a bottom quark correspondingly decrease. From Table 5.2 we can 
infer that xp decreases less quickly with the heavy meson mass than xs- Except for A^ 
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no empirical data for the masses of these baryons are known at present. These results 
for the mass differences among the bottom baryons are predictions of the model which 
can, in the future, be compared with experiment. As could have been inferred from 
the next to last row in Table 5.2 the absolute position of the bottom multiplet is about 
200 d= 50McV too low. On the absolute scale this apparently is only a 5% deviation 
from the data. Certainly a larger value o; ~ 1, which corresponds to a model for 
light vector resonance dominance of the heavy meson form factor [88], would yield 
an excellent agreement for the mass difference between Af, and the nucleon. On the 
other hand such a choice would slightly spoil the nice picture for the charm multiplet. 

The preceding calculations are based on the Nc —>■ oo limit in which the nucleon is 
infinitely heavy. Prom a common sense point of view this is peculiar since the nucleon 
is actually lighter than the heavy mesons being bound to it in the model. Hence, for 
comparison with experiment it is desirable to estimate kinematic effects associated 
with the nucleon's motion. These are expected [65] to lower the binding energy of the 
heavy baryons which have up to now come out too high (see Ac-N mass difference in 
Table 5.2, for example.). In order to estimate these kinematical effects in the bound 
state approach we have substituted the reduced masses 

into the bound state equations. In a non-relativistic treatment this corresponds to 
the elimination of the center of mass motion [65]. The results for the spectrum of 
the heavy baryons obtained from the replacement (5.44) are in displayed in Table 
5.4. Again we consider a as a free parameter. We notice that there is a remarkable 
improvement in the prediction for the Af, mass, which was previously the worst one. 
The changes in some of the mass parameters can approximately be compensated by 
a suitable re-adjustment of a. For a ~ 0.0 to —0.4 the agreement with the existing 
data is quite reasonable. When using the reduced meson masses the Ec baryon is 
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a 


0.0 


-0.1 


-0.2 


-0.3 


-0.4 


-0.5 


-0.6 


-0.7 


Expt. 


Up 


450 


469 


488 


508 


527 


546 


566 


585 




Xp 


0.212 


0.232 


0.246 


0.260 


0.273 


0.286 


0.299 


0.312 




UJs 


123 


134 


146 


158 


171 


184 


197 


210 




Xs 


0.410 


0.399 


0.387 


0.374 


0.361 


0.346 


0.331 


0.315 






158 


154 


151 


148 


145 


143 


140 


138 


168 


K 


221 


223 


225 


226 


227 


229 


230 


231 


233 


K 


342 


346 


353 


359 


363 


367 


371 


375 


308 


K 


460 


468 


475 


484 


490 


497 


505 


512 


? 


K* 


583 


587 


591 


596 


599 


601 


605 


607 


? 


N 


-1356 


-1338 


-1320 


-1302 


-1283 


-1265 


-1246 


-1228 


-1345 


Ab 


3285 


3282 


3280 


3278 


3275 


3272 


3271 


3269 


3356 ± 50 



Table 5.4: Parameters for heavy baryons and mass differences with re- 
spect to Ac. Primes indicate negative parity states, i.e. S-wave bound 
states. All energies are in MeV. In this calculation the reduced masses 
(5.44) enter the bound state equations from which the binding energies 
are extracted. The physical meson masses 1865MeV and 5279MeV are 
used when computing the mass differences to the nucleon and the A^ from 
these binding energies. Radially excited states are omitted because they 
are only very loosely bound, if at all. The empirical data are taken from 
the PDG [37], see also [92]. 
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a 


0.0 


-0.1 


-0.2 


-0.3 


-0.4 


-0.5 


-0.6 


-0.7 


Expt. 


UJp 


80 


94 


109 


124 


140 


155 


171 


188 




Xp 


0.346 


0.371 


0.394 


0.417 


0.439 


0.460 


0.479 


0.498 




E 


131 


126 


121 


117 


112 


108 


104 


100 


77 


E* 


235 


237 


239 


242 


244 


246 


248 


250 


269 


N 


-366 


-354 


-341 


-327 


-313 


-300 


-285 


-269 


-177 



Table 5.5: Same as Table 5.4 for even parity baryons in the kaon sector. 



always predicted a bit too light while it is too heavy when the physical masses are 
substituted in the bound state equation. For A'^ the situation is opposite. While 
the use of the physical meson masses gives too small a mass, the substitution of the 
reduced masses gives too large a prediction for the mass of this baryon. These results 
indicate that kinematical corrections are indeed important. 

It is interesting to see how far the heavy quark approach can be pushed to lighter 
quarks. To answer this question we have considered the strange quark. In the corre- 
sponding kaon sector the P-wave is only very loosely bound when the physical masses 
are substituted. On the other hand sizable binding energies arc obtained when the 
reduced masses are used [66]. This behavior is somewhat different from the charm 
and bottom sector and can be understood by noting that the difference M* — M is 
considerably reduced when using (5.44). In the heavy sectors (charm and bottom) 
this difference is small in any event. The resulting spectrum for the strange baryons is 
shown in Table 5.5. The comparison with the experimental data shows that even the 
use of the reduced masses does not provide sufficient binding. In the S-wave channel 
the situation is worse, even when the reduced masses are substituted bound states 
are not detected unless a < —1.0. The failure of the present approach in the strange 
sector strongly suggests that for these baryons a chirally invariant set-up [81] is more 
appropriate. 



Chapter 6 

Generalization of the Bound State 
Model 

6.1 Introduction 

Here we shall investigate the spectrum of excited states in the bound statc-soliton 
framework. Some aspects of this problem have already been treated [65,67,71,72,90]. 
We will deal with an aspect which does not seem to have been previously discussed in 
the literature. This emerges when one compares the excited heavy baryon spectrum 
with that expected in the constituent quark model (CQM) [73]. We do not have in 
mind specific dynamical treatments of the CQM but rather just its general geometric 
structure. Namely we shall just refer to the counting of states which follows from 
considering the baryon as a three body system obeying Fermi-Dirac statistics. We 
shall restrict our attention to the physical states for Nc = 3. In this framework the 
CQM counting of the heavy excited baryon multiplcts has been recently discussed [74]. 
At the level of two light flavors there are expected to be seven negative parity first 
excited A-type heavy baryons and seven negative parity first excited E-type heavy 
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baryons. On the other hand a similar counting [65,90] in the bound state treatments 
mentioned above yields only two of the A-type and five of the E-type. Thus there are 
seven missing first excited states. One thought is that these missing states should be 
unbound and thus represent new dynamical information with respect to the simple 
geometrical picture. There is certainly not enough data for the charmed baryons to 
decide this issue. However for the strange baryons there are ten established particles 
for these fourteen states. Hence it is reasonable to believe that these states exist for 
the heavy baryons too. In the CQM one may have two different sources of orbital 
angular momentum excitation; for example the relative angular momentum of the two 
light quarks, Lj and the angular momentum, Le of the diquark system with respect 
to the heavy quark. The parity of the heavy baryon is given by P = (— l)^^"*"^^. 
However, in the bound state models considered up to now there is only room for one 
relative angular momentum, r associated with the wave function of the heavy meson 
with respect to the soliton. The parity is given by P = (— 1)*^. Both models agree 
on the counting of the "ground" states (L/ — Le — r — 0). Also the counting of the 
states with (Lj — 0, Le — 1) agrees with those of r = 1 in the bound state model. 
However, the bound state model has no analog of the (L/ = 1, = 0) states and, 
in general, no analog of the higher Lj ^ states either. 

It is clear that we must find a way of incorporating a new angular momentum quan- 
tum number in the bound state picture. One might imagine a number of different 
ways to accomplish this goal. Here we will investigate a method which approximates 
a three body problem by an effective two body problem. Specifically we will consider 
binding excited heavy mesons with orbital angular momentum £ to the soliton. The 
excited heavy mesons may be interpreted as bound states of the original heavy me- 
son and a surrounding light meson cloud. Then the baryon parity comes out to be 
(— 1)'^"^^. This suggests a correspondence (but not an identity) r <-> L^, £ <-> L/ and 
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additional new states. An interesting conceptual point of the model is that it displays 
a correspondence between the excited heavy mesons and the excited heavy baryons. 

Almost immediately one sees that the model is considerably more complicated 
than the previous one in which the single heavy field multiplet H is bound to the 
soliton. Now, for each value of £ 7^ 0, there will be two different higher spin heavy 
multiplets which can contribute. In fact there is also a mixing between multiplets 
with different i, which is therefore not actually a good quantum number for the model 
(unless the mixing is neglected). 

Thus we will make a number of approximations which seem reasonable for an 
initial analysis. For one thing we shall neglect the light vector mesons even though 
we know they may be important. Wc shall also neglect the possible effects of higher 
spin light mesons, which one might otherwise consider natural when higher spin heavy 
mesons are being included. Since there is a proliferation of interaction terms among 
the light and heavy mesons we shall limit ourselves to those with the minimum number 
of derivatives. Finally, 1/M and nucleon recoil corrections will be neglected. The 
resulting model is the analog of the initial one used previously. Even though the 
true picture is likely to be more involved than our simplified model, we feel that the 
general scheme presented here will provide a useful guide for further work. 

We would like to stress that this bound state model goes beyond the kinematical 
enumeration of states and contains dynamical information. Specifically, the question 
of which states are bound depends on the magnitudes and signs of the coupling 
constants. There is a choice of coupling constants yielding a natural pattern of bound 
states which includes the missing ones. It turns out that it is easier to obtain the 
precise missing state pattern for the A-type heavy particles. Generally, there seem to 
be more than just the missing E-type heavy baryons present. However we show that 
the collective quantization, which is anyway required in the bound state approach. 
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leads to a splitting which may favor the missing heavy spin multiplets. 

This chapter is organized in the following way. Section 6.2 starts with a review of 
the CQM geometrical counting of excited heavy baryon multiplets. It continues with a 
quick summary of the treatment of heavy baryons in the existing bound state models. 
The comparison of the mass spectrum in the two different approaches reveals that 
there is a large family of "missing" excited states. This is discussed in general terms 
in section 6.3 where a proposal for solving the problem by considering the binding 
of heavy excited mesons to the Skyrmion is made. A correspondence between the 
angular momentum variables of the CQM and of the new model is set up. A detailed 
treatment of the proposed model for the case of the first excited heavy baryons is given 
in section 6.4. This includes discussion of the heavy meson bound state wave function, 
the classical potential energy as well as the energy corrections due to quantization of 
the collective variables of the model. It is pointed out that there is a possible way of 
choosing the coupling constants so as to bind all the missing states. The generalization 
to the excited heavy baryon states of arbitrary spin is given in section 6.5. This section 
also contains some new material on the interactions of the heavy meson multiplets 
with light chiral fields. Section 6.6 contains a discussion of the present status of 
the model introduced here. Finally, some details of the calculations are given in 
Appendices B.l and B.2. 



6.2 The Enigma of the Excited Missing States 

In this section, for the reader's convenience, we will briefiy discuss which heavy baryon 
states are predicted by the CQM as well as some relevant material needed for the 
bound state approach to the heavy baryon states. 

It is generally agreed that the geometrical structure of the CQM provides a reason- 
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able guide for, at least, counting and labeling the physical strong interaction ground 
states. When radial excitations or dynamical aspects are considered the model pre- 
dictions are presumably less reliable. In the CQM the heavy baryons consist of two 
light quarks (g) and a heavy quark (Q) in a color singlet state. Since the color singlet 
states are antisymmetric on interchange of the color labels of any two quarks, the 
overall wave function must, according to Fermi-Dirac statistics, be fully symmetric 
on interchange of flavor, spin and spatial indices. Here we will consider the case of 
two hght flavors. For counting the states we may choose coordinates [74] so that the 
total angular momentum of the heavy baryon, J is decomposed as 

J^Li + Le + S + Sh , (6.1) 

where Lj represents the relative orbital angular momentum of the two light quarks, 
Le the orbital angular momentum of the light diquark center of mass with respect 
to the heavy quark, S the total spin of the diquarks and Sh the spin of the heavy 
quark. In the "heavy" limit where the heavy quark becomes inflnitely massive Sh 
completely decouples. The parity of the heavy baryon is given by 

P5 = (-l)^^+^^ . (6.2) 

Since we are treating only the light degrees of freedom as identical particles it is only 
necessary to symmetrize the diquark product wave function with respect to the Lj, 
S and isospin / labels. Note that the diquark isospin I equals the baryon isospin. 
There are four possible ways to build an overall wave function symmetric with respect 
to these three labels: 



a) 


/ 


= 0, 




0, 


Li 


— even 


b) 


I 


= 1 , 




1 , 


Li 


= even 


c) 


I 


= 0, 


s = 


1 , 


Li 


= odd , 


d) 


I 


= 1 , 




0, 


Li 


= odd . 



(6.3) 
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There is no kinematic restriction on Le- 

Let us count the possible baryon states. The Lj — Le — heavy baryon ground 
state consists of Aq ( = i^) from a) and the heavy spin multiplet |Sq (|^) , Q 
from b). It is especially interesting to consider the first orbitally excited states. 
These all have negative parity with either {Le — 1, Lj — 0) or {Le = 0, Lj — 1). 
For Le — 1, a) provides the heavy spin multiplet |AqQ ^ , AqQ j| and b) 
provides EqQ-), {Sq (f ) , {Sq (f) , SgQ-)}. For = 1 c) 

provides Aq (i"), {Aq (i") , Aq (|~)}, {Aq (§") , Aq (§")}, while d) provides 
|Sq Q j , Sq (^I Altogether there are fourteen different isotopic spin multi- 
plets at the first excited level. The higher excited levels can be easily enumerated in 
the same way. For convenient reference these are listed in Table 6.1. 

It is natural to wonder whether all of these states should actually exist exper- 
imentally. This is clearly a premature question for the c and b baryons. However 
an indication for the first excited states can be gotten from the ordinary hyperons 
(or s baryons). In this case there are six well estabhshed candidates [75] for the 
A's [A(1405), A(1520), A(1670), A(1690), A(1800) and A(1830)]; only one |" state 
has not yet been observed. For the S's there arc four well established candidates 
[E(1670), E(1750), E(1775) and E(1940)]; two i" states and one |" state have not 
yet been observed. Thus it seems plausible to expect that all fourteen of the first 
excited negative parity heavy baryons do indeed exist. We might also expect higher 
excited states to exist. 

What is the situation in the bound state approach? To study this we shall briefly 
summarize the usual approach [65,67,90] to the excited heavy baryons in the bound 
state picture. In this model the heavy baryon is considered to be a heavy meson 

^We are adopting a convention where bold-faced angular momentum quantities are vectors and 
the regular quantities stand for their eigenvalues. 
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^« (i") 

{^Q(r).^Q(r)} 



= 1 



{AQ(r),A,(r)} 

{Ag(r),Acj(r)} 

{^^ (r).SQ(r)} 



L/ = 2n- 1 




2n - 




Table 6.1: Examples of the heavy baryon multiplets predicted by the 
CQM. 



112 



CHAPTER 6. GENERALIZATION OF THE BOUND STATE MODEL 



bound, via its interactions with the hght mesons, to a nucleon treated as a Skyrme 
sohton. The total baryon angular momentum is the sum (see Eq. (4.29)) 

J = g, + J-l + 5heavy , (6-4) 

where Sheavy is the spin of the heavy quark within the heavy meson. 

Now we can list the bound states of this model. First consider the r = states. 
According to Eq. (4.23), they have positive parity. Since Eq. (4.21) shows that 
A; = for binding, Eq. (4.19) tells us that the light "grand spin" g — 0. Equa- 
tion (4.29) indicates (noting I — J^°^) that there will be a Aq {^'^^ state as well as a 
|Sq (I^) , Sq (I^)} heavy spin multiplet. Actually the model also predicts a whole 
tower of states with increasing isospin. Next there will be an J = 2 heavy spin multi- 
plet with spins and parity | and | , and so forth. Clearly the isospin zero and one 
states correspond exactly to the Lj — Le — ground states of the constituent quark 
model. The isotopic spin two states would also be present if we were to consider 
the ground state heavy baryons in a constituent quark model with number of colors, 
Nc = 5. This is consistent with the picture [77] of the Skyrme model as a description 
of the large Nc limit. 

Next, consider the r — 1 states. These all have negative parity and (since the 
bound states have k = 0) light grand spin, g = 1. The J^°^ = 1 = choice yields 
a heavy multiplet |Aq ^ , Aq Q ^ | while the J^°^ = 1=1 choice yields the 
three heavy multiplets{EQ(i-)}, {Sq (f) , Eq (f)} and {Eq (f) ,Eq(|-)}. 
These three multiplets are associated with the intermediate sums \g + J^°'| = 0, 1, 2, 
respectively. It is evident that the seven states obtained have the same quantum 
numbers as the seven constituent quark states with L/ = and = 1. Proceeding 
in the same way, it is easy to see that the bound states with general r agree with 
those states in the constituent quark model which have Lj = and Le — r. This 
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may be understood by rewriting Eqs. (4.29) and (4.19) as 

J = r + + S-heavy , (6.5) 

where A: = for the bound states was used. Comparing this with the Lj = 
Umit of the constituent quark model relation (6.1) shows that there seems to be a 
correspondence 

'S'heavy SI h , 

r ^ Le , 

J'°' ^ S. (6.6) 

This correspondence is reinforced when we notice that / = J^°' in the bound state 
model and, for the relevant cases a) and b) in Eq. (6.3) of the constituent quark 
model, I — S also. We stress that Eq. (6.6) is a correspondence rather than an 
exact identification of the same dynamical variables in different models. It should 
be remarked that in the exact heavy and large Nc limits the heavy baryons for all 
values of r = g will have the same mass. When finite 1/M corrections are taken into 
account, there will always be, in addition to other things, a "centrifugal term" in the 
effective potential of the form g{g + l)/(2M|a;p), which makes the states with larger 
values of g, heavier. It should also be remarked that the above described ordering 
of heavy baryon states in the bound state approach applies only to the heavy limit, 
where Shcavy decouples. For finite heavy quark masses, multiplets are characterized 
by the total grand spin g + (Sheavy Then states like Aq j and Aq ^| j no longer 
constitute a degenerate multiplet. 
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The Planetary Conjecture 



It is clear that the bound state model discussed above contains only half of the 
fourteen negative parity, first excited states predicted by the CQM. The states with 
Lj are all missing. Since the enumeration of states in the CQM was purely 
kinematical one might at first think that the bound state model (noting that the 
dynamical condition k = was used) is providing a welcome constraint on the large 
number of expected states. However, experiment indicates that this is not likely to be 
the case. As pointed out in the last section, there are at present good experimental 
candidates for ten out of the fourteen negative parity, first excited ordinary hyperons. 
Thus the missing excited states appear to be a serious problem for the bound state 
model. 

The goal of the present work is to find a suitable extension of the bound state 
model which gives the same spectrum as the CQM. Reference to Eq. (6.1) suggests 
that we introduce a new degree of freedom which is related in some way to the light 
diquark relative angular momentum Lj. To gain some perspective, and because we are 
working in a Skyrme model overall framework, it is worthwhile to consider the heavy 
baryons in a hypothetical world with Nc quark colors. In such a case there would be 
Nc — 1 relative angular momentum variables and we would require Nc — 2 additional 
degrees of freedom. Very schematically we might imagine, as in Fig. 6.1, one heavy 
meson H and Nc — 2 light mesons A4i orbiting around the nucleon. One might 
imagine a number of different schemes for treating the inevitably complicated bound 
state dynamics of such a system. Even in the Nc — S case it is much simpler if we can 
manage to reduce the three body problem to an effective two body problem. This 
can be achieved, as schematically indicated in Fig. 6.2, if we link the two "orbiting" 
mesons together in a state which carries internal angular momentum. The "linked 
mesons" will be described mathematically by a single excited heavy meson multiplet 
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Figure 6.1: Schematic planetary picture for large Nc excited heavy 
baryons in the bound state approach. 
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Figure 6.2: Schematic picture of the "two body" approximation for the 
Nc — 3 excited heavy baryons. 

field. One may alternatively consider these "linked mesons" as bare heavy mesons 
surrounded by a light meson cloud. Such fields are usually classified by the value^, i of 
the relative orbital angular momentum of a qQ pair which describes it in the CQM. 
We will not attempt to explain the binding of these two mesons but shall simply 
incorporate the "experimental" higher spin meson fields into our chiral Lagrangian. 
Different i excitations will correspond to the use of different meson field multiplets. 
Prom now on we will restrict our attention to A^c 3. 

Taking the new degree of freedom £ into account requires us to modify the previous 
formulas describing the heavy baryon. Now the parity formula (4.23) is modified to 

Pb = (-1)'+'' , (6.7) 

which is seen to be compatible with the CQM relation (6.2). Now Eq. (4.29) holds 

^Actually if we want to picture the linked mesons as literally composed of a meson-meson pair, 
we should assign relative orbital angular momentum ^ — 1 to these bosonic constituents and allow 
for both light pseudoscalars and vectors. 
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but with the hght grand spin g modified to, 

g = r + K' . (6.8) 

Note that K in Eq. (4.18) has been incorporated in 

K' — /light + <S^iight + ^ ■ (6.9) 

The new correspondence between the bound-state picture variables and those of the 
CQM is: 



heavy 


<— > 


Sh 


r 


<— > 


Le 


£ 


<— > 


Li 


, jsol 


<— > 


S . 



lught + ^ught + ^ S. (6.10) 

Previously /ught + 'S'ught = K had zero quantum numbers on the bound states; now 
the picture is a little more complicated. We will see that the dynamics may lead to 
new bound states which are in correspondence with the CQM. Equation (6.10) should 
be interpreted in the sense of this correspondence. 

It is easiest to see that the lowest new states generated agree with the CQM for 
£ = even, which corresponds to negative parity heavy mesons. In this case k = 
or equivalently k' — i may be favored dynamically. Then the last line in Eq. (6.10) 
indicates that J^°\ which can take on the values and 1, corresponds to the hght 
diquark spin S in the CQM. This leads to the CQM states of type a) and b) in 
Eq. (6.3). This is just a generalization of the discussion for the ground state given 
in section 6.2. Now let us discuss how the states corresponding to c) and d) can be 
constructed in the bound state scenario. Apparently we require i — odd, i. e. positive 
parity heavy mesons. For 7 = we also have J^°' = 0. Hence the last fine in Eq. (6.10) 
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requires A; = 1 for 5* = 1. To generate states of type d) also k = 1 would be needed 
in order to accommodate / = J^°^ — 1 and S — 0. Actually for the case A; = 1 and 
— 1 states with -S" = 0, 1, 2 would be possible. The states with = 1, 2 should be 
ruled out by the dynamics of the model. 

One may perhaps wonder whether we are pushing the bound state picture too far; 
since things seem to be getting more complicated why not just use the constituent 
quark model? Apart from the intrinsic interest of the soliton approach there are 
two more or less practical reasons for pursuing the approach. The first is that the 
parameters of the underlying chiral Lagrangian are, unlike parameters such as the 
constituent quark masses and inter-quark potentials of the CQM, physical ones and 
in principle subject to direct experimental test. The second reason is that the bound 
state approach actually models the expected behavior of a confining theory; namely, 
when sufficient energy is applied to "stretch" the heavy baryon it does not come 
apart into a heavy quark and two light quarks but rather into a nucleon and a heavy 
meson. The light quark-antiquark pair which one usually imagines popping out of 
the vacuum when the color singlet state has been suitably stretched, was there all the 
time, waiting to play a role, in the bound state picture. The model may therefore be 
useful in treating reactions of this sort. 



6.4 Model for the Missing First Excited States 

Before going on to the general orbital excited states it may be helpful to see how the 

dynamics could work out for explaining the missing seven Aq and type, negative 

parity, excited states. In the new bound state picture these correspond to the choices^ 

^Actually, i was introduced for convenience in making a comparison with the constituent quark 
model. It is really hidden in the heavy mesons which, strictly speaking, are specified by the light 
cloud angular momentum Jught and parity. We can perform the calculation without mentioning i. 
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£ = 1, r — 0. As discussed, we are considering that the orbital angular momentum 
i is "locked-up" in suitable excited heavy mesons. As in Eq. (4.17), r appears as 
a parameter in the new heavy meson wave-function. The treatment of the excited 
heavy mesons in the effective theory context, has been given already by Falk and 
Luke [78]. For a review see [79]. The case (for orbital angular momentum=l) where 
the light cloud spin of the heavy meson is 1/2 is described by the heavy multiplet 

^ - ^ ^ {S + 757.^'^) , (6.11) 

where S is the fluctuation field for a scalar (J^ = 0+) particle and satisfying 
Vij,A^ — 0, similarly corresponds to an axial {J^ — 1+) particle. The case where the 
light cloud spin is 3/2 is described by 



dfJ-iy c^lv il/J, ^) 



(6.12) 



satisfying the Rarita-Schwinger constraints H/j.'f'^ — H/jV^ — 0. The field Tn^, — 7]^^ 
(with Vi'Ti,^ = T;;* = 0) is a spin 2 tensor (J^ = 2+) and (with Vi^Bi" = 0) is 
another axial (J^ = 1"*"). Currently, experimental candidates exist for the tensor and 
an axial. 

In order to prevent the calculation from becoming too complicated we will adopt 
the approximation of leaving out the light vector mesons. This is a common approx- 
imation used by workers in the field but it should be kept in mind that the effect of 
the light vectors is expected to be substantial. 

The kinetic terms of the effective chiral Lagrangian (analogous to the first term 
of Eq. (4.11)) are: 

^Ckin = -iMV^'Tr lUDf^H] -iMV^'Tr lU'^D^Ha] , (6.13) 



where is a characteristic heavy mass scale for the excited mesons. For simphcity^ 
*A more general approach is to replace A4 on the right hand side of Eq. (6.13) by the same M used 
in Eq. (4.11) and to add the splitting terms -2M{Ms - M)Tr [HH] - 2M{Mt - M)Tr [W^W**] . 
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we are neglecting mass differences between the £ = 1 heavy mesons. The interaction 
terms involving only the Tl and fields, to lowest order in derivatives, are 



(6.14) 



These generalize the second term in Eq. (4.11) and ds, d^ and /st (which may be 
complex) are the heavy meson-pion coupling constants. Similar terms which involve 
£ 7^ 1 multiplets are not needed for our present purpose but will be discussed in the 
next section. 

As in section 6.2, the wave-functions for the excited heavy mesons bound to the 
background Skyrmion are conveniently presented in the rest frame V = 0. The 
analogs of Eq. (4.14) become 



ftk 





(6.15) 



and ('^^o) ~^ 0. Now the wave-functions in Eq. (6.15) are expanded as: 

u{\x\) 



fih — 

fa — 
J i,lh 



u{\x\) 



(6.16) 



where u stands for a sharply peaked radial wave-function which may differ for the 
two cases. Other notations are as in Eq. (4.15). Note that the constraint = 

implies that 

{a,\,^,,i,d^O. (6.17) 

It is interesting to see explicitly how the extra angular momentum £ = 1 is "locked- 
up" in the heavy meson wave-functions. For the Ti. wave-function, the fact that 
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•flight — + 'S'light takes the value 1/2 leads, using Eq. (6.9), to the possible values 
/c' = or 1. The corresponding wave-functions are 

<^>M{k' = k', = 0) = ^ , $,,(fc' = fc3 = l) = ^, (6.18) 



where, for the present case, we are taking r — 0. For the Hi wave-function it is 
important to satisfy ji = | Jughtl =3/2 condition (6.17). This may be accomplished 
by combining with suitable Clebsch-Gordan coefficients an i = 1 wave-function with 
the S'light = 1/2 spinor to give 

^i,id{k' = K = 2) = wl^'^6n6di , 

^i,id {k' = A;^ = 1) = ^w\-^'^5n5a2 - - l=wP5a5ai , (6.19) 

where wf^^^ = [Sji ± iSj2) and wf'^ = is a spherical decomposition. 

The main question is: Which of the channels contain bound states? Note that, 
for the reduced space in which x ■ r has been removed as in Eq. (6.16), k' is a 
good quantum number. Furthermore, because the wave-function M(|a;|) is sharply 
peaked, the relevant matrix elements are actually independent of the orbital angular 
momentum r. The classical potential for each k' channel may be calculated by setting 
r — and substituting the appropriate reduced wave-functions from Eqs. (6.18) and 
(6.19) into the interaction Lagrangian (6.14). (see Appendix B.l for more details.) 
The A;' = channel gets a contribution only from the ds term in Eq. (6.14) while the 
k' — 2 channel receives a contribution only from the term. On the other hand, all 
three terms contribute to the k' — 1 channel. The resulting potentials are: 

Vik' = 0) = -IdsF'iO), (6.20) 
V(k'^2) = -IdrF'iO), (6.21) 

y (.' = !)= I ™ (W) Lf I^^(O). (6.22) 

{n,\v\n) {n,\v\n,) j \ ^ji^ I^t 
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The classical criterion for a channel to contain a bound state is that its potential 
be negative. Since -F'(O) > we require for bound states in the A;' = and k' — 2 
channels 

ds>0 , dT>0 , (6.23) 

respectively^. For bound states in the k' — 1 channel we must examine the signs 
of the eigenvalues of Eq. (6.22). Assuming that Eq. (6.23) holds (as will be seen to 
be desirable) it is easy to sec that there is, at most, one k' = 1 bound state. The 
condition for this bound state to exist is 



(6.24) 



The (primed) states which diagonalize Eq. (6.22) are simply related to the original 
ones by 




cos 9 sin 9 
y — p* sin 9 p cos 9 J 

tan2.= ^^l^- 



(6.25) 



(6.26) 



where p is the phase of /st- ^ and $i are shorthand notations^ for the appropri- 
ate wave-functions. Clearly, the results for which states are bound depend on the 
numerical values and signs of the coupling constants. At the moment there is no 
purely experimental information on these quantities. However, it is very interesting 
to observe that if Eqs. (6.23) and (6.24) hold, then the missing first excited Aq states 
are bound. To see this note that the heavy baryon spin is given by Eq. (4.29) with g 



Hn a more general picture where i = 3 excited heavy mesons are included, the k' = 2 channel will 

also be described by a potential matrix. Then the criterion for is modified. (See next section.) 
^Strictly speaking, to put on a parallel footing to ^i,id we should replace 

^id \/i (-P^^^)ife;;;' ('^fc)dd'*;'d' witli tlie spin 3/2 projection operator, (-P^/^)^^.;;, = 

I [SikSw - ^ejik (see Appendix B.l). 
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defined in Eqs. (6.8) and (6.9). For the Ag-type states, noting that / = J^°^ = in 
the Skyrme approach gives the baryon spin as 

J = g + Sheavy (Ag states) . (6.27) 

The r — choice enables us to set g — k'. With just the three attractive channels 
A;' = 0, A;' = 1 and k' — 2 we thus end up with the missing first three excited Aq 
heavy multiplets Aq (i") , {Aq (|") , Aq (§")} and {Aq (§") , Aq (|~)}. It should 
be stressed that this counting involves dynamics rather than pure kinematics. For ex- 
ample, it may be seen from Eqs. (6.20)-(6.22) that it is dynamically impossible to have 
four bound heavy multiplets {k' — 0, k' — 2 and two k' — 1 channels). The missing 
first excited Eg-type states comprise the single heavy multiplet |Sq Q ^ , Sq ^ |. 
At the classical level there are apparently more bound multiplets present. However, 
we will now see that the introduction of collective coordinates, as is anyway required 
in the Skyrme model [80] to generate states with good isospin quantum number, will 
split the heavy multiplets from each other. Thus, deciding which states are bound 
actually requires a more detailed analysis. 

We need to extend Eq. (4.24) in order to allow the i = 1 heavy meson fields to 
depend on the collective rotation variable A{t): 

n{x, t) = A{t)n^{x) , n^{x, t) = A{t)n^,{x) , (6.28) 

where He and Hie are given in Eq. (6.15). Note, again, that the matrix A{t) acts on the 
isospin indices. We also have Hoc — due to the rest frame constraint V'^H/j^ — 0. 
Now substituting Eq. (6.28) as well as the first of Eq. (4.24) into the heavy field 
Lagrangian^ yields [59] the collective Lagrangian^ 

Le,n^la''n''-x{k')K'-n , (6.29) 
''Note that Eq. (6.13) contributes but Eq. (6.14) does not contribute. 

^In Eq. (6.29) k' is defined to operate on the heavy particle wave-functions rather than on their 
conjugates. This is required when the heavy meson is coupled to the Skyrme background field since 
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where O is defined in Eq. (4.25) and is the Skyrme model moment of inertia. In the 
vector meson model the induced fields (po and cui) are determined from a variational 
approach to a^. The quantities x {k') are given by (see Appendix B.2). 



X{k') 



k' ^0 

i(3cos2^-l) , (6.30) 



where the angle 9 is defined in Eq. (6.26). (Note that if hght vector mesons are 
included the expressions for x would be more involved as the induced fields will also 
contribute.) In writing Eq. (6.30) it was assumed that the first state in Eq. (6.25) 
(i.e. $' rather than $'j) is the bound one; the collective Lagrangian is constructed as 
an expansion around the bound state solutions. We next determine from Eq. (4.28), 
the canonical (angular) momentum J^°^ as a'^ Q — x {k') K'. The usual Legendre 
transform then leads to the collective Hamiltonian 

H,on^^{r''' + xik')K'f . (6.31) 

Again we remark that J^°^ = /. It is useful to define the light part of the total heavy 
baryon spin as 

j = r + K' + J^°^ , (6.32) 

and rewrite Eq. (6.31) as 

^^coii = ^ [(1 - X (^0) I' + X {k') {j -rf + x {k') ix {k') - 1) K''] . (6.33) 

The mass splittings within each given k' multiplet due to i^cou are displayed in Ta- 
ble 6.2. This table also shows the splitting of the k' multiplets from each other due to 
Aq is made as {qqq) (qQ) rather than {qqq) (Qq). For convenience in Eqs. (4.15) and (6.16) we have 
considered the conjugate wave-functions (since they are usual in the hght sector). This has been 
compensated by the minus sign in the second term of Eq. (6.29). 



6.4. MODEL FOR THE MISSING FIRST EXCITED STATES 



125 



k' 



K' + J 



sol 



V 



a'H, 



coll 



Candidates for r = 

missing states 




1 
2 



-Ids F'iO) 



A 

idTF'iO) 



A, 



Q (l 



X 

_3_ 
16 



{AQ(i-),AQ(r)} 

{AQ(r)>AQ(i-)} 



1 


1 
2 
1 
2 
3 



-IdsF'iO) 
A 



(x 
(x- 



idTF'iO) 



-if + x 

l)' + 3x 

16 



N(r).s«(r)}, 
{E«(r),E«(r)}3 



15 
16 

27 
16 



Table 6.2: Contributions to energies of new predicted 
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ative binding potential in the k' — 1 channel. Furthermore x — x(l) 
Eq. (6.30); it satisfies ^ 



2- 



126 



CHAPTER 6. GENERALIZATION OF THE BOUND STATE MODEL 



the classical potential in Eqs. (6.20)-(6.22). Note that the slope of the Skyrme profile 
function F'(0) is of order 1 GeV. The coupling constants ds, d^, /st, based on d ~ 0.5 
for the ground state heavy meson, are expected to be of the order unity. Hence the 
binding potentials V are expected to be of the rough order of 500 MeV. The inverse 
moment of inertia is of the order of 200 MeV which (together with —\ <X^\) 
sets the scale for the "1/A^c" corrections due to i^cou- As mentioned before, if the 
couphng constants satisfy the inequalities (6.23) and (6.24), all the Aq multiplets 
shown will be bound. At first glance we might expect all the Eq states fisted also to 
be bound. However the H^ow corrections increase as / increases, which is a possible 
indication that many of the Eq's might be only weakly bound. In a more complete 
model they may become unbound. Hence it is interesting to ask which of the three 
displayed candidates for the single missing Eg multiplet is mostly tightly bound in 
the present model. Neglecting the effect of V we can see that iJcoU raises the energy 
of candidate 3 less than those of candidates 1 and 2. Furthermore, for the large range 
ofx, —\ <x< 1 — candidate 3 suffers the least unbinding due to i?coii of any of 
the I — 1 heavy baryons listed. The Aq states suffer still less unbinding due to i^cou- 



6.5 Extension to the Higher Orbital Excitations 

We have already explicitly seen that the "missing" first orbitally excited heavy baryon 
states in the bound state picture might be generated if the model is extended to also 
include binding the first orbitally excited heavy mesons in the background field of a 
Skyrme soliton. Prom the correspondence (6.10) and associated discussion we expect 
that any of the higher excited heavy baryons of the CQM might be similarly generated 
by binding the appropriately excited heavy mesons. In this section we will show in 
detail how this result can be achieved in the general case. An extra complication. 
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which was neglected for simphcity in the last section, is the possibility of baryon 
states constructed by binding heavy mesons of different £, mixing with each other. 
For example {r — 1 , i — 0} type states can mix with {r — 1 , £ — 2} type states, other 
quantum numbers being the same. Since r + £ must add to 1, this channel could not 
mix with {r = 1 , £ = A}. An identical type of mixing - between {Le = 1 , Lj = 0} 
and {Le — I , Lj — 2} - may also exist in the CQM. The present model, however, 
provides a simple way to study this kind of mixing as a perturbation. 

To start the analysis it may be helpful to refer to Table 6.3, which shows our 
notations for the excited heavy meson multiplct "fluctuation" flclds. The straight 
i^'s contain negative parity mesons and the curly 7^'s contain positive parity mesons. 
Further details are given in Ref . [78] . Note that each field is symmetric in all Lorentz 
indices and obeys the constraints 



(6.34) 



as well as for 'H^^...ij_^. The general chiral invariant interaction with the lowest number 
of derivatives is 

£d + A + >Cg , (6.35) 

where 



The final piece. 



n=0 



n=0 



A = M^/p,Tr fi/^--'^»p'^75^^,..^ 



n=0 



+ h.c. 



n=0 



+ h.C. . 



£g = M ^ ^„Tr \W-^-T^l,l,H,,...A + h.c. 



n=0 



(6.36) 



(6.37) 



128 



CHAPTER 6. GENERALIZATION OF THE BOUND STATE MODEL 



field 



Ji 



J' 



H 



1/2 



0-, 1- 



n 



1/2 
3/2 



0+, 1+ 
1+, 2+ 



H„ 



H, 



2 
2 



3/2 
5/2 



1- , 2- 

2- , 3- 



^ = even £-1/2 (£-l)-,r 

^w-M. ^ = even i + 1/2 r , (£ + 1)" 

£ = odd £-1/2 (£-l)+,£+ 

£ = odd £+1/2 £+, (£+1)+ 



Table 6.3: Notation for the heavy meson multiplets. ji is the angular 
momentum of the "light cloud" surrounding the heavy quark while is 
the spin parity of each heavy meson in the multiplet. 
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exists in general, but does not contribute for our ansatz. Terms of the form 



Tr 



H 



Tr 



n 



(6.38) 



can be shown to vanish by the heavy spin symmetry In the notation of Eq. (6.14), 
<^s = d^Oi c?T = c^si ^-iid /sT = /so- A new type of coupUng present in Eq. (6.36) also 
connects multiplets to others differing by A£ = ±2. These are the terms with odd 
(even) n for H (7i^)-type fields. The interactions in Eq. (6.37) connecting multiplets 
differing by A£ = ±1 turn out not to contribute in our model. In the interest of 
simplicity we will consider all heavy mesons to have the same mass. This is clearly 
an approximation which may be improved in the future. 

The rest frame ansdtze for the bound state wave functions which generalize Eq. (6.15) 
are (note ji — n + 1/2): 



h---i„,lh 



ii---in,lh 




Jl 



1 

2 ' 



(6.39) 



with identical structures for ^ — > 7Y. Note that again a, I, h represent respectively 
the isospin, light spin and heavy spin bivalent indices. Extracting a factor of A • r as 
we did before in Eqs. (4.15) and (6.16) leads to 



"'ii---i„,lh 



U\\X\ 



(6.40) 



with similar notations. The relevant wave-functions are the ipij^...i„,di {k', k'^, r). k' was 
defined in Eq. (6.9); we will see that it remains a good quantum number. Since the 
terms which connect the positive parity {H type) and negative parity (7i type) heavy 
mesons (Eq. (6.37)) vanish when the ansdtze (6.39) are substituted, the baryon states 
associated with each type do not mix with each other in our model. We thus list 



130 



CHAPTER 6. GENERALIZATION OF THE BOUND STATE MODEL 



separately the potentials for each type. For the £ = even baryons (associated with H 



mesons) , 



V [k' = 0] 

V [k' 7^ 0] 



-2^poi^'(0), 



F(0) 



k' dp(^k'~i) 



-1 

2 



-/p(fc'-i) 



-(-1) 



fc/ 2k' + 3 dpk' 



2k' + 1 2 



(6.41) 



while for the £ — odd baryons (associated with H mesons), 



V [k' = 0] 

V [k' 7^ 0] 



-2^soi^'(0) , 



F'(0) 



(-1) 



k' ds{k'^i) 



(-1) 



g/s(fc'-l) 

u 2k' + 3 



Sfc' 



(6.42) 



2A;' + 1 2 

Details of the derivations of Eqs. (6.41) and (6.42) are given in Appendix B.l. The 
ordering of matrix elements in Eqs. (6.41) and (6.42), for a given k', is such that the 
first heavy meson has a light spin, ji = k' — ^ while the second has ji = k' + \. The 
H type {H. type) channels with k' — even (odd) involve two mesons with the same 
£ — k'. The H type {Ti. type) channels with k' — odd (even) involve two mesons 
differing by A£ = 2, i.e., £ — k' — 1 and £ — k' + 1. This pattern is, for convenience, 
illustrated in Table 6.4. Also shown, for each k', are the number of channels which 
are expected to be bound according to the CQM. 

It is important to note that Table 6.4 holds for any value of the angular momentum 
r, which is a good quantum number in our model. For the reader's orientation, we 
now locate the previously considered cases in Table 6.4. The standard "ground state" 
heavy baryons discussed in section 6.2 are made from the H meson with i = and 
ji — 1/2. They have r — and k' — 0. The seven negative parity heavy baryons 
discussed in section 6.2 also are made from the H meson with £ — and ji — 1/2. 
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k' 31 


if mesons 
^ # 


7i mesons 
^ # 


1/2 


1 


1 1 


3/2 






2 


1 

1 

1 


3/2 
2 ' 

5/2 


2 

1 

2 


1 

2 

3 


5/2 

3 

7/2 


2 



4 


3 

1 

3 



Table 6.4: Pattern of states for Eqs. (6.41) and (6.42). Note that ji = n+\ 
is the hght cloud spin of the heavy meson. The columns marked ^ stand 
for the number of channels which are expected to be bound, for that 
particular k' , according to the CQM. 
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They still have k' — 0, but now r = 1. The seven "missing" first excited heavy 
baryons discussed in section 6.4 have r — and are made from the i — 1, Ti. and 
Hn mesons with ji — 1/2 and ji — 3/2. There should appear one bound state for 
fc' = 0, one bound state for k' = 1 and one bound state for fc' = 2 in the "?Y-meson" 
section of Tabic 6.4. Note that the number of states expected in the CQM model for 
k' — 2 is listed in Table 6.4 as two, rather than one. In the absence of Ai — 2 terms 
connecting H/j, and H^u (see the last term in Eq. (6.36)) i would be conserved for 
our model and only the i — 1 state would be relevant. This was the approximation 
we made, for simplicity, in section 6.4. The other entry would have i = 3 and would 
decouple. When the Ai = 2 mixing terms are turned on, the £ = 1 and i = 3, k' = 2 
channels will mix. One diagonal linear combination should be counted against the 
Lj — 1 CQM states and one against the Lj — 3 CQM states. 

To summarize: for the H-tjpe mesons, the even k' channels should each have 
one bound state, while the odd k' channels should have none. The situation is very 
different for the 7i-type mesons; then the even k' ^ channels should contain two 
bound states while the odd k' channels should contain one bound state. The A;' = 
channel should have one bound state. 

For the iJ-type meson case, the pattern of bound states mentioned above would 
be achieved dynamically if the coupling constants satisfied: 



These follow from requiring only one negative eigenvalue of Eq. (6.41) for k' = even 
and none for k' — odd. Similarly requiring for the 7i-type meson case in Eq. (6.42), a 
negative eigenvalue for k' — 0, one negative eigenvalue for k' — odd and two negative 



rfpo > , 





{k' = odd) . 



(6.43) 



6.5. EXTENSION TO THE HIGHER ORBITAL EXCITATIONS 



133 



eigenvalues for A;' > and even leads to the criteria, 



dso > 



(-1) 



k' 



ds{k'-i)dsk' 



'2k' + 3\ _ 8 
2k' + l) ~ 3 



/s(ifc'-l) 



> , {k' >0) 



f 2k' + 3\ 

ds{k'-i) + I Trrrri- 1 ^Sk' > , {k' ^ even ^ 0) . 



(6.44) 



2k' + 1^ 

From Eqs. (6.43) and (6.44) it can be seen that all the ci's are required to be positive. 
Furthermore these equations imply that the |/|'s which connect heavy mesons with 
Ai — 2 are relatively small (compared to the d's) while the |/|'s which connect heavy 



mesons with A£ — are relatively large. In detail this means that /p 



'(fc'-i) 



should be 



small for odd k' and large for even k' with just the reverse for fs{k'-i) ■ This result 
seems physically reasonable. 

As in the example in the preceding section we should introduce the collective 
variable A{t) in order to define states of good isospin and angular momentum. This 



again yields some splitting of the different 



K' + J 



sol 



members of each k' bound state. 



Now, each k' channel (except for k' = 0) is described by a 2 x 2 matrix. Thus there 
will be an appropriate mixing angle 9, analogous to the one introduced in Eq. (6.25), 
for each k' and parity choice {i.e., H-type or H-type field). The collective Lagrangian 
is still given by Eq. (6.29) but, in the general case. 



X±{k') 



1 



^±(^k' + ^]cos2e 



2k'{k' + 1) 

In this formula the different signs correspond to the two possible eigenvalues, 

,k'-i / 



(6.45) 



d(k'-i) - 



2k' 



2k' + 1 



dy + 



32 „ 2 
/(fe'-i) 



F'(0) 



(6.46) 



of the potential matrix. For example, referring to Table 6.4, we would expect the 
k' — 2, Ti-type meson case to provide two distinct bound states and hence both 
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X+(2,H) and x_(2,?i) would be non-zero. On the other hand, we would expect no 
bound states in the k' — 3, H-type meson case so x±(3, H) should be interpreted as 
zero. 

It is convenient to summarize the energies of the predicted states in tabular form, 
generalizing the example presented in Table 6.2. The situation for baryons with 
parity = —(—1)'' ("^^-type mesons) is presented in Table 6.5. For definiteness we 
have made the assumption that the constraints (6.44) above are satisfied. In order 



= J^°l 


k' 


K' + J^°i 


V 




Candidates for r = 
missing states 





In - 1 


2n- 1 


A+ 


n(2n — l)xi 


{A((2n-3/2)-) , A((2n-l/2)-)} 


In 


2n 


A+ 

A_ 


n(2n + l)x^ 

/K2» + l),Yl 


{A((2n-l/2)-) , A((2n + l/2)-)} 


1 


2n- 1 


2n-2 
2n- 1 
2n 


A+ 


n{2n - l)x+ + 1 - 2nx+ 

n(2n-l)x+ + l-X+ 
n(2n - \)x\ + 1 + (2n - l)x+ 


{S((2n-3/2)-) ,S((2n-l/2)-)}i 


2n 


2n- 1 

2n 
2n + 1 


A+ 


n(2n + \)x\ + 1 - (2n + l)x+ 
n(2n + \)x\ + 1 - X+ 
n(2n + l)x+ + 1 + 2nx+ 


{S((2n-3/2)-) , S((2n-l/2)-)}2 
{S((2n + l/2)-),S((2n + 3/2)-)}3 


2n- 1 

2n 
2n + 1 




n(2n + l)x- + 1 - (2n + l)x- 
n(2n + l)x^ + 1 - X- 
n(2n + \)^_ + 1 + 2nx- 


{S((2n-3/2)-),E((2n-l/2)-)}4 

{S((2n+l/2)-) ,S((2n + 3/2)-)}5 



Table 6.5: Contributions to energies of the new predicted states made 
from 7Y-type heavy mesons. Note that n is a positive integer. The n = 
case is given in Table 6.2. The A+ entries in the V column are more 



tightly bound than the A_ entries. 



K' \J 



sol 



is the light part of the 



heavy baryon angular momentum for r = (See Eq. (6.32).). 
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to explain Table 6.5 let us ask which states correspond to the {Lj = 3, Le — 0) 
states in the CQM. Reference to Table 6.1 shows that three negative parity A-type 
heavy multiplets and one negative parity E-type heavy multiplet should be present. 
The correspondence in Eq. (6.10) instructs us to set r = and, noting Eq. (6.9) , to 
identify 

K' + J^°i ^Li + S . (6.47) 

The A-type particles are of type c) in Eq. (6.3) so we must take S — 1. Hence, since 
jsoi _ Q fQj. ^_type particles, we learn that k' can take on the values 2, 3 and 4. For 
k' — 2, the second line of the k' column yields two possible multiplets (energies A+ 
and A_) with n = 1 and structure |A Q ^ , A Q Wc should choose one of these 
to be associated with (L/ = 3, = 0) and the other with {Lj = 1, Le = 0) in 
the CQM. We remind the reader that £ is not a good quantum number so that the 
correspondence £ ^ Lj in Eq. (6.10) only holds when the Ai — 2 mixing terms are 
neglected. For fc' = 3, the first hue of the k' column correctly yields one multiplet 
with n = 2 and structure {A (^| j , A For k' = 4, the second line of the k' 

column yields two multiplets with n — 2 and structure |A Q j , A One of 

these is to be associated with (Lj — 3, Le — 0) and the other with (L/ — 5, Le — 0) 
in the CQM. Now let us go on to the E-type heavy multiplets. These are of type d) in 
Eq. (6.3) and yield S = 0. Hence K'+ J'°' ^ Lj and K' + J'°' = 3. Five candidates 
for this {e(| ) , S(| )} multiplet are shown in the last column of Table 6.5. These 



K' + J 



sol 



column. 



consecutively correspond to the choices n = 2,2, 1,2, lin the 
As before it is necessary for an exact correspondence with the CQM that one of these 
should be dynamically favored (much more tightly bound) over the others. Again, 



note that the choice 



K' + J 



sol 



3 does not uniquely constrain the value of i. 
Next, the situation for baryons with parity = (— 1)*^ (-f^~type baryons) is presented 
in Table 6.6.. For definiteness we have made the assumption that the constraints 
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I 

= J^°l 


k' 


K' + J^°' 


V 


X Hcou 


Candidates for r = 
missing states 





2n 


2n 


A+ 


n{2n + l)x\ 


{A((2n-l/2)+) , A((2n+l/2)+)} 


1 


2n 


2n-l 

2n 
2n + 1 


A+ 


n{2n-l)xl + l-{2n+l)x+ 
n{2n + l)x\ + 1-X+ 
n{2n + l)x+ + 1 + 2nx+ 


{S((2n-3/2n,S((2n-l/2)+)}, 
{S((2n-l/2)+) , S((2n + l/2)+)}2 
{S((2n + l/2)+),S((2n + 3/2)+)}3 



Table 6.6: Contributions to energies of the new predicted states made 
from if-type heavy mesons. Other details as for Table 6.5. 



(6.43) above are satisfied. This eliminates the odd k' states and agrees with the CQM 
counting. For example, we ask which states correspond to the (L/ = 2, Le — 0) 
states in the CQM. Reference to Table 6.1 shows that one positive parity A-type 
heavy multiplet and three positive parity E-type heavy multiplets should be present. 
For r = we have the correspondence K' + J^°^ ^ Lj + S. The A-type particles are 
of type a) in Eq. (6.3) so we must set k' = 2. The first line in Table 6.6 then yields, 
with n — 1 the desired |A {^^^ , A [^^^ } heavy multiplet. The E particles are of 
type b) in Eq. (6.3) so that ^K' + J^°^| can take on the values 1, 2 and 3. The last 
three lines in Table 6.6, with n = 1, give the desired multiplets: |e (^^^^ , E (|^)}, 
|e (^I^^ , E (f^)} and |e (I"*") , E (i^)}- In this case all the states should be 
bound so that the splittings due to i^cou are desired to be relatively small. The 
present structure is simpler than the one shown in Table 6.5 for the 7:^-type cases. 
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6.6 Conclusions for the Generalized Heavy Baryon 
Model 

In this chapter [91] we have pointed out the problem of getting, in the framework of 
a bound state picture, the excited states which are expected on geometrical grounds 
from the constituent quark model. We treated the heavy baryons and made use of 
the Isgur-Wisc heavy spin symmetry. The approach may also provide some insight 
into the understanding of light excited baryons. The key problem to be solved is 
the introduction of an additional "source" of angular momentum in the model. It 
was noted that this might be achieved in a simple way by postulating that excited 
heavy mesons, which have "locked-in" angular momentum, are bound in the back- 
ground Skyrmion field. The model was seen to naturally have the correct kinematical 
structure in order to provide the excited states which were missing in earlier models. 

An important aspect of this work is the investigation of which states in the model 
are actually bound. This is a complicated issue since there are many interaction 
terms present with a priori unknown coupling constants. Hence, for the purpose of 
our initial investigation we included only terms with the minimal interactions of the 
hght pseudoscalar mesons. The large M limit was also assumed and nucleon recoil 
as well as mass splittings among the heavy excited meson multiplets were neglected. 
We expect, based on previous work, that the most important improvement of the 
present calculation would be to include the interactions of the light vector mesons. 
It is natural to expect that possible interactions of the light higher spin mesons also 
play a role. In the calculation of the ground state heavy baryons the light vectors 
were actually slightly more important than the light pseudoscalars and reinforced the 
binding due to the latter. Another complicating factor is the presence, expected from 
phenomenology, of radially excited mesons along with orbitally excited ones. 
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It is interesting to estimate which of the first excited states, discussed in sec- 
tion 6.4, are bound. The criteria for actually obtaining the missing states in the model 
with only light pseudoscalars present are given in Eqs. (6.23) and (6.24). Based on 
the use of chiral symmetry for relating the coupling constants to axial matrix ele- 
ments and using a quark model argument to estimate the axial matrix elements, Falk 
and Luke [78] presented the estimates (their Eqs. (2.23) and (2.24)) — 3cis — d 
and |/st| = -^d. With these estimates Eqs. (6.23) and (6.24) are satisfied. Note 
that d > provides binding for the ground state heavy baryons. However we have 
checked this and find that, although we arc in agreement for |/st| we obtain instead 
d^ = 3ds = —d. Assuming that this is the case then it is easy to see that the only 
bound multiplet will have k' — 1. This leads to the desired E-type multiplet and 
one of the three desired A-type multiplets being bound, but not the A;' = and 2, 
A-type multiplets. Clearly, it is important to make a more detailed calculation of 
the light meson-excited heavy meson coupling constants. We also plan to investigate 
the effects of including light vector mesons in the present model. It is hoped that the 
study of these questions will lead to a better understanding of the dynamics of the 
excited heavy particles. 

Finally we would hke to add a few remarks on studies of the excited "light" 
hyperons within the bound state approach to the SU(3) Skyrme model. In that model 
the heavy spin symmetry is not maintained since the vector counterpart of the kaon, 
the K*, is omitted; while the kaons themselves couple to the pions as prescribed by 
chiral symmetry. On the other hand the higher orbital angular momentum channels 
{i.e. r > 2) have been extensively studied. The first study was performed by the 
SLAG group [82]. However, they were mostly interested in the amplitudes for kaon- 
nucleon scattering and for simplicity omitted flavor symmetry breaking terms in the 
effective Lagrangian. Hence they did not find any bound states, except for zero modes. 
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These symmetry breaking terms were, however, included in the scattering analysis 
of all higher orbital angular momentum channels by Scoccola [83]. The only bound 
states he observed were those for P- and S-waves. After collective quantization these 
are associated with the ordinary hypcrons and the A(1405). As a matter of fact these 
states were ah^cady found in the original study by Callan and Klcbanov [59]. It is clear 
that the orbital excitations found in the bound state approach to the Skyrme model 
should be identified as the £ — states. Furthermore when the dynamical couphng 
of the collective coordinates (^4, f2) is included in the scattering analysis [84] the only 
resonances which are observed obey the selection rule \ J —l/2\ < r < |J+1/2|, 
where r denotes the kaon orbital angular momentum. This rule is consistent with 
£ — in our model. In order to find states with £ ^ in this model one would also 
have to include pion fluctuations besides the kaon fluctuations for the projectile-state. 
As indicated in the previous sections, these fluctuating fields should be coupled to 
carry the good quantum number i. The full calculation would not only require this 
complicated coupling but also an expansion of the Lagrangian up to fourth order in the 
meson fluctuations off the background soliton. Such a calculation seems impractical, 
indicating that something like our present approximation, which treats these coupled 
states as elementary particles, is needed. 



Appendix A 
Part I Appendix 



A.l Scattering kinematics 

The general partial wave scattering matrix for the multichannel case can be written 
as: 

Sab^5ab + 2iTal, . (A.l) 

For simplicity, the diagonal isospin and angular momentum labels have not been 
indicated. 

By requiring the unitarity condition S^S — 1 one deduces for the two-channel 
case the following relations: 

ImiTn) = \Tn\' + \T2i\\ 

Im{T22) = |^22p + |Tl2|^ (A.2) 
Im{Tu) = T*,Tu + T*^T22 , 

where T12 — T21. In the present case we will identify 1 as the tttt channel and 2 as 
the KK channel. In order to get the relations between the relative phase shifts and 

140 



A.l. SCATTERING KINEMATICS 



141 



the amplitude we need to consider the following parameterization of the scattering 
amplitude: 

, ^ , (A.3) 

where S^^k — ^tv + and < 77 < 1 is the elasticity parameter. By comparing 
eq. (A.3) and eq. (A.l) one can easily deduce: 

7^2 = l-4|ri2|' . (A.4) 

Analogously, for Taa we have: 

where / and / label the angular momentum and isospin respectively. Extracting the 
real and imaginary parts via 



T3l ^ 



2 



J l-Tjl cos(2(5^.,) ^ ^ ^ 

C;i = ' 3 ^ (A.6) 

leads to the very important bounds 

\Ria;l\<l, 0</i;,<l. (A.7) 

Unitarity also requires 1^/2.^1 < 1/2 . 

Now wc relate these partial wave amplitudes to the invariant amplitudes. The 
invariant amplitude for iTiipi) + T^jip2) —>■ T^kips) + ''^lip^) is decomposed as: 

5ij5kiA{s, t, u) + 8ik8jiA{t, s, u) + 8ii8jkA{u, t, s) , (A. 8) 

where s, t and u are the usual Mandelstam variables. Note that the phase of eq. (A. 8) 
corresponds to simply taking the matrix element of the Lagrangian density of a four- 
point contact interaction. Projecting out amplitudes of definite isospin yields: 

r{\(s, t, u) = 3A(s, t, u) + A{t, s, u) + A(u, t, s) , 
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Tl,{s,t,u) 
Tl{s,t,u) 



A{t, s, u) — A{u, t, s) , 
A{t,s,u) +A{u,t,s) . 



(A.9) 



The needed / = tttt — > KK amplitude can be obtained as: 



TUs,t,u) = -V6Ai7i'{p,)7i%p2),K+{p3)K-{p,)) . 



(A.IO) 



We then define the partial wave isospin amplitudes according to the following formula: 
Tab-A^) = IVP^b l\cosePi{cos9)T^,{s,t,u) , (A.ll) 



where 9 is the scattering angle and 



Pa = 



' s — Ami 



S 16% 



- 0{s - Ami) ■ 



(A.12) 



is a symmetry factor which is 2 for identical particles (tttt case) and 1 for distin- 
guishable particles [KK case) . 



A. 2 Unregularized amplitudes 

A. 2.1 Amplitudes for the tttt tttt channel 

The current algebra contribution to A{s,t,u) is 

(s-ml) 



Aca{s,t,u) = 2- 



F3 



The amplitude for the vectors can be expressed in the following form 



t{u — s) u(t — s) 
mi — t mi — u 



where Qpt^t^ is the coupling of the vector to two pions. 



(A.13) 



(A.14) 
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For the scalar particle we deduce 



9x2 



Af^{s,t,u) 



7o js - 2nr^) 
2 rnj^ — s 



(A.15) 



To calculate the tensor exchange diagram we need the spin 2 propagator [44] 

(A.16) 



g2 - mj^ 



where 



e 



-9nu + 



m 



/2 



(A.17) 



A straightforward computation then yields the /2 contribution to the tttt scattering 
amplitude: 



Af^{s,t,u) = 



7i 



2(m 



/2 



16 4 10 2 1 2 I 1 /i2 I 2\ 

— —ml + —m^s - -s + - {t + u ) 



2 nri^s^ 

3 6m^2 6^/2, 



+ 



(A.18) 



A.2.2 Amplitudes for tt^tt'^ K+K 



Current algebra amplitude: 



2F2 ■ 

TT 



Vector meson contribution: 



Av^ector(7r°7r°,X+X-) = 



9K*Kn 



K 

2 ™2\2 



t{s ~ u) ^ u{s — t) 
m\, — t m\;* — u 



^ " XTrix* - 1 m\, - 



1 1 

+ 



Direct-channel contribution for the scalar: 

1 



(s-2ml){s-2ml) 



m%-s 



(A.19) 



(A.20) 



(A.21) 
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Cross-channel contribution for the scalar: 

{ml + iii'l - t)- {ml + ml - u)^ 



AKs{7r\',K+K-) 



^-2 

I KT.Ktt 



8 



m_fj-* t 



m 



u 



.(A.22) 



Direct channel tensor contribution: 



2{m% - s) 



4mi ^ 2 



t (ml + mlY 



+ 



u {ml-\-m\y 



4m% + 2 



3 \4m% 2 
Cross-channel tensor contribution: 



2- - o + "^TT 



(A.23) 



AK*(7r'7r',K+K-) = 
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{t-ml- ml) - — ^(t^ - {ml - mlY) 
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+ {t < — > u) . 
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Appendix B 
Part II Appendix 



B.l Classical Potential 



Here we will show how to compute the relevant matrix elements associated with the 
classical potential. 

For any fixed value of /c' 7^ the heavy meson light cloud spin (Jught) takes 
the values ji — k' ^ ^ since K' — Jught + -^light , where /ught is the heavy meson 
isospin. Hence the classical potential will be, in general, a 2 x 2 matrix schematically 
represented as 



V {k' 0) 



'-1) i-^fil-fJ'k' l^l-^Ml-'-Mfc' 



(B.l) 



Here corresponds to the ji = k' — \ state while corresponds to 

ji — k' + \. In order to compute the potential there is no need to distinguish even 
parity heavy mesons Ti. from odd parity ones H. The diagonal matrix elements are 
obtained by substituting the appropriate rest frame ansatz (6.39) into the general 
potential term as: 
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= 4 ^ (-1)" / dn ri,...i„,ai (^'' ^3, r) a,v ■ Taa'^,,...,^4'i' (k\ r) , (B.2) 

where ji = n + and n = k' ^1 for the two diagonal matrix elements. The operator 
which mesures the total light cloud spin ji is 

("flight),- , , = ® ^iih <^i„i„ + ® i-i^ahn) ® 5i2i2 ® ■ ■ ■ ® ^wn 

+ • • • + • • • ® (-ie„i„^-J . (B.3) 

where Caij is the totally antisymmetric tensor. The isospin operator is 

i'light = ^ . (B.4) 

We can write Eq. (B.3) compactly in the following way 

flight = s + i , (B.5) 

where s = y- Due to the total symmetrization of the vectorial indices we have I — n. 
We want to stress that s and I do not necessarily agree with Sought and Indeed 
for associated with V. in Eq. (6.16), I — and Jught — s — Sught + (■ while for 
associated ^i^id with I — 1. Now we have, for fixed n — ji — \^ the following 
useful result: 

■ / ■ ^ IX / dn-il) JiightV' = 77- 

By using Eq. (B.6) we can write Eq. (B.2) as 

{-ly-dn — : — / d^lip* {k',k3,r) J light- 1 hght'ip{k',k3,' 

= (-l)"c^„^[fc'(fc' + l)-i,(j, + l)-^" . 

For ji = /c' =p I we get the diagonal matrix elements for both, the H type as well as 
the H type fields 

F'(V\\ I dk'-i , 3i = k' — \ , 

(-if-^-l , /2.' + 3X . ^. (B.8) 



j ^.^^ ^ jdnr{s-Jii^u)i^ f fdn V* JiightV' ■ (B.6) 

J 7/.7Z + 1 J 2.7/7 



(B.7) 
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where we used n = ji — 1/2. 

For the non-diagonal matrix elements we consider the contribution to the potential 
due to the following / type term: 



2 

This corresponds to the transition between ji — n+ ^ and ji — n+ ^ states. Now we 
notice that by construction any wave function ip must satisfy the condition 

where P^/^ is the spin 3/2 projection operator 

The condition (B.IO) yields the following identity 

/ dQi^l.,^^,, ik',K,r)Ti,, (P^/^)^.^^^^>,,..,„MT ik',K,r) . (B.12) 
Using the fact that P^^^r commutes with K', we get 

i^'^')jk;ii' ^dd'i'n-ir.A'v {k\ k',, t) = Ar^i,...„,,,z {k\ k',, r) , (B.13) 
where A?" is a normalization constant. It is evaluated as 

W = / {k\k',,r)T:,,, (P='/2)^^^^^,r,V^n...„,<^'T {k^k'.^r) = ^ . (B.14) 
The non-diagonal matrix element is, up to a phase factor in Eq. (B.9) 

^/nF'(0)y|, VfcVO. (B.15) 

For A;' = we have only one diagonal element with ji — |. The second line of 
Eq. (B.8) provides 

V{k' = 0) = ~F'{0)do . (B.16) 
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B.2 Collective Lagrangian 

Here the relevant matrix elements associated with the collective coordinate Lagrangian 
are computed. We will restrict k' to be nonzero since there is no contribution for k' — 
to the collective Lagrangian. 

The kinetic Lagrangian for H type and H type fields is: 



n n ~ 

(B.17) 

In the following we will not distinguish between the H and Ti. types of field. We 
need to consider the collective coordinate Lagrangian for a given k' classical bound 
channel in the heavy meson rest frame. For k' the bound state wave-function 
can schematically be represented as 

\Bound State; k') = a \H^,...^^,J + (5 \H^,...^^,) , (B.18) 

where jap + = 1. 

The collective coordinate Lagrangian (SLcoii), induced by the heavy meson kinetic 
term, is obtained by generalizing Eqs. (4.24) and (6.28) to the higher excited heavy 
meson fields, introducing the collective coordinate A{t) rotation via 

H,,...,^{x,t) = A{t)H,,...U^) , (B.19) 

where the Hi-^...i^c{x) classical ansatz is given in Eq. (6.39). The contribution for fixed 
k' j^O is: 



+ \P\'J dni^l..,^,^,,{k',k'„r) 'fA...,,„,'iik',k'„r) 
= -lal'Jdnr {k\ k'^, 31 ^k'- 1/2) a . JughtV' {k', k'^, ji ^ k' - 1/2) 
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J dnr {k', K,3i + 1/2) n ■ /light {k', K,3i ^k' + 1/2) , 

(B.20) 

where the over all minus sign in Eq. (B.20) is required, as explained in section 6.4. 
According to the Wigner-Eckart theorem: 



k'{k' + l)-Mji + l) + 



j dni/j^K'i/j , (B.21) 



2k'{k' + l) 

we thus obtain the following heavy meson contribution to the collective coordinate 

Lagrangian for k' ^ 

SL^n = -X{k') n K' . (B.22) 

The quantity x(/c') is given by 

y(k') = — -i \- ± (k' + -\ cos 2^1 , (B.23) 

' 2k'{k'^l) [2 V 27 J ' ^ ' 

where |q;|^ — = ±cos2^ was used. In Eq. (B.23) the ± sign corresponds to the 
two possible eigenvalues in the potential matrix for given k' ^ 0. 
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